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ABSTRACT 


report  deals  with  the  following  topics: 

1)  Detection  in  the  presence  of  interference.  Earlier  studies  of 
analogue  detection  jn  an  interference  dominated  environment 
were  reported  in  Volune  III.  Here  the  emphasis  is  on  detectors 
using  digital  techniques.  Costs  of  clipping  and  sampling  are 
examined  for  optimal  detectors  as  veil  as  for  more  readily 
realizable  suboptimal  instrumentations.  Some  attention  is  given 
to  the  problem  of  recovering  part  of  the  clipping  loss  thro'^h 
multilevel  quantization. 

2)  Adaptation  to  certain  noise  field  properties.  The  report  examines 
detectors  which  adjust  the  detection  threshold  in  accordance  with 
various  measured  noise  parameters.  An  instrumentation  using 
measured  noise  power  in  setting  the  threshold  possesses  non- 
parametric  properties  that  make  the  fclse  alarm  rate  asymptotically 
independent  of  the  noise  amplitude  distribution.  A  second  (digital) 
instrumentation  achieves  reduced  dependence  on  noise  spectral 
properties  by  adjusting  the  threshold  in  accordance  with  the  count 
of  zeros  at  each  hydrophone. 

3)  Active  Sonar.  One  possible  detector  for  active  sonar  signals 
consists  of  a  pair  of  widely  separated  receivers  whose  outputs 
are  cross-correlated.  The  performance  of  such  a  receiver  is 
analyzed  for  noise  environments  dominated  by  reverberation  or  by 
ambient  noise.  The  effect  of  linear  frequency  modulation  in  the 
transmitted  signal  is  examined. 

A)  Vertical  directionality  of  ambient  noise.  Numerical  computations 
are  performed  to  estimate  '■he  directionality  of  ambient  noise  at 
shallow  depths  for  various  realistic  velocity  profiles.  The  basic 
noise  model  is  that  proposed  by  Talham. 
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FOREWORD 


This  report  is  the  fourth  in  a  scries  describing  work  performed  by  Yale  University 
under  subcontract  to  Electric  Boat  division  of  General  Dynamics  Corporation.  The 
report  covers  the  period  l  July  1965  to  i  July  1906.  An  unclassified  supplement  to 
this  volume  has  been  bound  as  a  separate  document,  Report  No.  U417-07-084. 
Electric  Boat  is  prime  contrrctor  of  the  SUBIC  (Submarine  Integrated  Control) 
Program  under  Office  of  Na^al  Research  contract  NOnr  2512(01)).  UCDR.  E.W.  Lull, 
USN,  is  Project  Officer  for  ONR;  J.  W.  Herring  is  Project  Manager  for  Electric 
Boat  division  under  t.ie  direction  of  Dr.  A.  J.  van  Woerkom. 
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I . _ Introduction 

The  following  is  a  summary  of  work  performed  under  contract 
8050-31-55001  between  Yale  University  and  the  Electric  Boat  division  during 
the  period  1  July  1965  to  1  July  1966.  More  detailed  discussions  of  the 
results  as  well  as  their  derivations  are  contained  in  a  series  of  six 
progress  reports  that  are  appended.  Four  of  these  represent  continuations 
of  investigations  on  passive  detection  reported  in  earlier  volumes  of  this 
series,  one  initiates  a  new  effort  in  the  direction  of  active  sonar  signal 
processing  (continued  in  subsequent  reports)  while  the  final  one  is 
concerned  with  a  special  topic  in  noise  propagation. 

One  of  the  central  themes  in  Volume  III  of  this  series  was  detection 
in  a  noise  field  dominated  by  a  plane  wave  interference.  The  performance 
of  the  optimal  (likelihood-ratio)  detector  was  compared  with  that  of  a 
conventional  detector  and  found  to  be  substantially  better  in  a  strong 
interference  environment.  Several  suboptimal  analogue  procedures  (null 
steering)  were  proposed  and  analyzed.  In  the  present  volume  this  effort 
is  continued  with  the  examination  of  procedures  suitable  for  digital 
implementation.  The  simple  introduction  of  sampling  and  hard  limiting  into 
the  previously  analyzed  suboptimal  instrumentations  was  found  to  lead  to 
degradations  In  performance  which  grow  with  interference  to  ambient  noise 
ratio.  A  study  was  therefore  undertaken,  a)  to  separate  the  loss  into  a 
sampling  component  and  a  clipping  component,  b)  to  determine  what  fraction 
of  the  loss  was  independent  of  the  particular  implementation  (hence  a  basic 
clipping  loss)  ,  and  c)  to  estimate  the  improvement  attainable  through  use 
of  n  rather  than  2  level  quantization. 

In  most  analytical  treatments  of  the  passive  detection  problem  the 
statistical  properties  of  signal  and  noise  are  assumed  to  be  known  a  priori. 
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This  is  clearly  unrealistic  in  most  practical  sonar  problems.  As  one 
attempts  to  relax  the  assumption  one  is  lead  Lo  instrumentations  which 
measure  relevant  statistical  parameters  and  adapt  the  detector  to  the 
measurements .  Depending  on  one's  point  of  view,  the  resulting  system 
might  be,  at  one  extreme,  a  configuration  adapting  to  the  complete  space 
time  structure  of  the  noise  field  or,  at  the  other  extreme,  a  system 
adapting  only  to  one  or  two  parameters  deemed  particularly'  important  by 
the  designer.  The  null  steering  procedures  mentioned  earlier  may  be 
said  to  fall  into  the  second  category,  the  single  unknown  parameter  in 
this  case  being  the  interference  bearing.  Another  parameter  almost 
certainly  unknown  a  priori  is  the  total  noise  power  or  signal-to-noise 
ratio.  Since  it  critically  affects  the  detection  threshold,  measurement 
and  adaptation  suggests  itself.  Initial  studies  concerning  thiB  problem 
were  contained  in  Volume  III.  Two  reports  in  the  present  volume  extend 
this  effort  by  developing  procedures  insensitive  not  only  to  noise  power 
level,  but  also  noise  amplitude  distribution  and,  to  some  extent,  noise 
spectral  properties. 

The  initial  effort  in  the  active  sonar  field  deals  with  a  detector 
cross-correlating  the  outputs  of  two  spatially  separated  receivers  in  an 
environment  that  may  be  dominated  by  either  ambient  noise  or  reverberation. 
One  of  the  primary  aims  of  this  study  was  the  development  of  on  analytical 
framework  for  the  discussion  of  reverberation  that  could  be  carried  over 
to  future  investigations  of  active  signal  design  and  processing. 

II.  Detection  in  the  Presence  of  Interference 

Two  reports  in  this  volume  deal  with  detection  in  a  noise  environment 
dominated  by  a  strong  plane  wave  interference.  Both  assume  that  the  output 
of  each  hydrophone  is  converted  into  digital  (generally  binary)  form  prior 
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to  further  processing.  Fi_port  No.  26  examines  a  digital  version  of  the 
n  ill  steering  procedure  analyzed  in  Volume  ITT  (Report  To.  21).  The 
hydrophone  outputs  are  sampled  and  hard-limited,  digitally  delayed  for 
alignment  with  the  interference,  subtracted  pairwise  for  interference 
elimination,  after  which  beamforming  on  the  target  and  detection  is 
accomplished  in  conventional  manner.  The  sampling  procedure  degrades 
performance  by  two  mechanisms. 

1)  If  the  sampling  period  corresponds  to  the  minimum  delay 
increment  available,  then  the  sampling  rate  fixes  the  number 
of  beams  that  can  be  formed.  Unless  the  interference  falls 
precisely  on  the  axis  of  one  of  these  beams  it  cannot  be 
eliminated  perfectly  by  the  null  steering  mechanism. 

2)  Even  if  the  interference  lies  precisely  on  a  beam  axis  (and 
could  therefore  be  eliminated  perfectly  by  an  analogue  nulling 
procedure),  hard  limiting  spreads  the  spectrum  over  an  infinite 
frequency  interval  (even  if  the  input  spectrum  is  bandlimited) . 
Therefore  any  finite  sampling  rate  leads  to  a  certain  amount  of 
spectral  foldovcr  with  attendant  loss  of  informat io*. 

The  first  effect  is  easily  studied  by  considering  r  sampled  analogue 

version  of  the  null  steering  detector.  Signal  interference  aud  ambient 

noise  are  assumed  prewhitened  to  a  bandwidth  w  tad/sec  so  that  the 

o 

Nyquist  rate  is  vQ/it  samples/sec.  The  residual  interference  power  is  a 
maximum  when  the  interference  is  located  midway  between  the  two  wojL 
closely  adjacent  beams.  In  this  worst  case  the  sampling  rate  f 
(samples/sec)  required  to  hold  the  residual  interference  to  a  power  no 
larger  than  that  of  the  ambient  noise  is  given  approximately  by 
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f 

e 


0.9 


(Nyquist  rntt) 


for 


I/N  >.  1 


(1) 


I/N  is  Che  interference  to  noise  ratio.  Thus,  for  an  interference  to 
noise  ratio  of  100  (20  db) ,  one  must  sample  at  about  9  times  Nyquist 

rate . 

Having  isolated  the  beam  steering  problem,  one  can  now  deal  with  the 
clipping  loss  and  the  additional  sampling  loss  due  to  mechanism  2)  by 
constraining  the  interference  to  lie  on  the  beam  axis.  In  order  to 
separate  the  effects  of  clipping  ard  sampling  as  much  as  possible  (always 
keeping  in  mind  that  there  would  be  no  type  2)  sampling  loss  in  the' 
absence  of  clipping)  it  is  convenient  to  consider  two  cases: 

a)  Sampling  at  the  Nyquist  rate  of  the  prewhitened  signal,  noise 
and  inference  processes,  so  that  the  samples  are  statistically 
independent . 

b)  Sampling  at  a  rate  faster  than  the  above  by  a  sufficient  amount 
so  that  the  output  -,ignal-to-noise  ratio  (considered  as  a 
function  of  sampling  rate)  closely  approaches  its  asymptotic 
value . 

In  case  a)  one  finds  that  the  clipping  and  sampling  loss  (ratio 
of  output  signal-to-noise  ratio  without  clipping  to  output  signal -to-noise 
ratio  with  clipping)  is  given  roughly  by 


clipping  and  sampling  loss  (in  db)  =*  10  log  1.4  ~  J  for  I/N>1  (2) 

The  above  result  is  essentially  independent  of  the  number  of  hydrophones 
once  the  array  size  exceeds  a  relatively  small  minimum.  The  improvement 
of  the  null  steering  detector  over  the  equivalent  (clipped  and  Nyquist 
rate  sampled)  conventional  detector  is  quite  substantial  for  strong 
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interferences.  For  a  AO  element  linear  array  at  1/N  »  100  (20  db) 
the  improvement  amounts  to  16.6  db  of  output  signal-to-noise  ratio. 

This  figure  increases  by  roughly  6  db  for  every  10  db  increment  in 
interference  to  noise  ratio,  but  is  relatively  insensitive  to  changes 
in  the  number  of  hydrophones. 

It  is  most  convenient  to  express  the  results  of  part  b)  (fast 
sampling)  in  terms  of  those  of  part  a)  (Nyquist  rate  sampling)  and  an 
improvement  factor  I.  The  latter  rises  linearly  wita  sampling  rate  until 
it  approaches  an  asymptotic  value  near  tht  sampling  rate  given  by  Eq .  (1)  . 
Thus  sampling  above  that  rate  brings  little  'benefit.  An  estimate  of  the 
true  cost  of  clipping  can  be  obtained  by  subtracting  from  Eq .  (2)  the 
asymptotic  value  of  I  (expressed  in  db) .  The  result  (for  reasonably 
large  arrays)  is  a  loss  of  about  2  db  for  I /N  ■*  1  and  an  increase  of 
roughly  2  db  for  every  10  db  increase  in  interference  to  noise  ratio. 
Similar  computations  for  n  set  of  non-white  spectra  indicate  that  these 
results  do  not  depend  critically  on  spectral  shape. 

Report  No.  26  evaluated  the  cost  of  clipping  for  a  particular 
instrumentation  (null  steering)  „  It  therefore  leaves  open  the  question 
whether  the  very  appreciable  losses  encountered  with  high  interference  to 
noise  ratios  are  peculiar  to  this  particular  system  or  whether  they  are 
a  fundamental  property  of  all  clipped  ins trumentatione.  To  ahed  light  on 
this  question,  Report  No.  28  considers  a  likelihood  ratio  detector 
operating  on  clipped  hydrophone  outputs.  The  results  indicate  that,  for 
large  interference  to  noise  ratios,  the  optimum  system  is  a  clipped 
analogue  of  Anderson's  null  ste  .ring  scheme  (beam-form  on  interference, 
subtract  average  value  of  this  interference  estimate  from  each  hydrophone, 
then  bean  form  on  the  target).  The  performance  difference  (large  I/M) 
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between  this  system  and  the  null  steering  scheme  of  Report  No.  26  is  minor 
(less  than  1  db  ) .  Thus  the  clipping  coat  determined  in  Report  No.  26 
appears  to  be  basic  and  not  subject  to  reduction  by  improved  processing 
procedures . 

An  interesting  physical  insight  into  the  reasons  for  the  clipping 
loss  can  be  gained  from  the  following  argument.  Consider  a  linear  ar  y 
mechanically  steering  on  target  so  that  the  signal  wave-front  reaches  all 
hydrophones  simultaneously.  The  situation  is  illustrated  in  the  diagram. 
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The  vert_eal  dimension  corresponds  to  spatial  displacement  along  the 
array,  the  horizontal  to  tine  displacement  (successive  samples)  of  the 
output  of  one  hydrophone.  In  the  absence  of  interference  the  noise  is 
independent  from  hydrophone  to  hydrophone.  The  signal  wavefront  causes 
one  polarity  to  be  preferred  along  each  vertical  line  and  it  is  this 
preponderance,  however  slight,  which  allows  detection  when  a  sufficient 
number  of  samples  (vertical  lines)  is  examined.  In  the  presence  of  a 
strong  interference  the  polarities  of  samples  along  an  interference 
wavefront  are  all  the  same  unless  the  interference  time  function  is  near 
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zero  at  th^  particular  instant  represented  by  that  wavefront.  A  set  of 

samples  whose  polarity  is  completely  determined  by  the  interference 

contributes  no  infor nation  concerning  the  presence  or  absence  of  the 

signal.  Detection  must  then  be  accomplished  on  the  basis  of  the  residual 

samples,  those  corresponding  to  sufficiently  low  instantaneous  values  of 

interference  so  that  their  polarity  is  affected  by  signal  and  ambienc 

noise.  It  is  a  simple  natter  to  calculate  the  probability  that  all 

samples  along  an  interference  wavefront  have  the  same  polarity.  For 

large  interference  to  noise  ratios  this  probability  differs  from  unity 

by  a  factor  of  the  order  of  \J N/I  .  Hence  only  about  ^/'N/I  of  the  total 

samples  are  useful  in  signal  detection.  More  quantitatively,  one  finds 

from  this  line  of  reasoning  that  the  clipping  cost  for  slow  sampling 

is  asymptotically  10  log  ^  J  [compare  with  Eq.  (2)]. 

in-  v  N 

Report  No.  28  also  examines  the  improvement  attainable  through  use 
of  multilevel  quantization  in  place  of  hard  limiting.  Approximate 
compute1. tions  indicate  that,  for  an  interference  to  noise  ratio  cf  20  db, 
the  use  of  4  quantization  levels  reduces  the  clipping  loss  from  12  db 
to  9  db  for  slow  sampling  and  from  6  db  to  db  for  fast  sampling. 

More  than  8  levels  are  required  to  reduce  the  fast  sampling  loss  to  3  db . 
One  must  note  also  that  the  adjustment  of  multilevel  quantizers  becomes 
more  critically  dependent  on  the  ambient  noise  power  ns  the  number  of 
levels  increases. 

III.  Adaptation  to  Certain  Noise  Field  Propci  ties 

s 

The  problem  of  detection  in  a  noise  environment  of  unknown  power 
level  was  first  considered  in  Volume  III  (Report  No.  18).  An  instrumenta¬ 
tion  was  there  proposed  which  measured  the  a’  erage  noise  power  over  the 
observation  interval  and  used  the  result  to  set  the  detector  threshold. 
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In  a  Gaussian  noise  environment  and  with  arrays  of  more  than  a  very  small 
number  of  hydrophones  the  resulting  system  performed  almost  as  well  as 
one  having  full  a  priori  knowledge  of  the  noise  power.  Report  No.  23  of 
the  present  series  examines  the  performance  of  the  same  adaptive  threshold 
detector  in  more  general  noise  environments.  One  interesting  conclusion 
is  that  the  false  alarm  rate  becomes  independent  of  the  noise  probability 
density  as  the  observation  time  increases.  In  this  very  important  sense, 
therefore,  the  detector  is  asymptotically  nonparamctric .  The  report  also 
examines  the  dependence  of  detection  probability  (fixed  false  alarm  rate) 
on  the  noise  amplitude  distribution.  With  Gaussian  noise  there  is  a 
degradation  cf  performance  equivalent  to  the  loss  of  about  *5  hydrophone 
(relative  to  the  perfor;,  .ncc  of  a  conventional  detector  operating  in  a 
noise  environment  of  known  power).  On  the  other  hand,  it  is  not  difficult 
to  construct  examples  of  non-Gaussian  noise  in  which  the  detection  proba¬ 
bility  (fixed  false  alarm  rate)  exceeds  that  of  a  conventional  detector. 
This  appears  to  be  true  especially  with  noise  of  an  impulsive  character. 

Report  No.  28  discusses  an  adaptive  threshold  for  a  clipped  conven¬ 
tional  detector  (PCA  detector) .  Clipping  eliminates  the  dependence  on 
noise  power,  but  there  remains  a  dependence  on  noise  spectral  properties. 
In  particular,  the  fact  that  changes  in  noise  bandwidth  affect  the 
dependence  between  adjacent  samples  suggests  that  improved  performance 
night  be  obtiined  by  adjusting  the  detector  threshold  in  accordance  with 
the  noise  bandwidth.  An  approximate  measure  of  this  bandwidth  can  be 
obtained  from  a  count  of  zeros  at  the  output  of  each  hydrophone.  The 
report  considers  the  performance  of  a  detector  which  ties  its  threshold 
to  the  zero  count.  The  functional  dependence  cf  the  threshold  on  the  zero 
count  is  based  on  a  nominal  spectrum  anc  the  performance  of  the  detector 
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is  then  investigated  for  a  wide  variety  of  different  noise  spectra.  The 
procedure  is  by  no  means  perfect,  but  the  false  alarm  rate  turns  out  to 
be  far  less  dependent  on  spectral  properties  than  that  of  a  fixed 
threshold  detector. 

IV.  Active  Sonar 

Report  No.  27  analyzes  an  active  sonar  detector  using  cross-correlation 
between  the  signals  received  at  two  spatially  separated  receivers.  The 
transmitted  signal  consists  of  a  single  pulse  of  linearly  frequency 
modulated  carrier.  The  pulse  duration  is  large  compared  with  the  carrier 
period.  Transmitter,  target,  receiver  and  all  scattering  centers  are 
assumed  to  be  stationary,  so  that  no  Doppler  shifts  occur.  The  received 
signal  is  assumed  to  be  a  delayed,  but  otherwise  perfect  ,  replica  of  the 
transmitted  signal  anu  a  similar  assumption  is  made  concerning  the  returns 
from  various  scattering  centers  which  make  up  the  reverberation.  This 
description  of  the  sonar  return  is  obviously  far  too  idealized  for  many 
purposes.  it  ignores  most  of  the  important  transformations  experienced 
by  the  signal  in  a  realistic  propagation  channel.  However,  it  has  the 
merit  of  not  obscuring  the  effects  of  the  reverberation  model  whoso 
development  and  study  is  one  of  the  primary  purposes  of  the  report. 

The  reverberation  is  attributed  to  the  effect  of  a  large  number  of 
omnidirectional  point  scatterers,  Toisson  distributed  throughout  the 
illuminated  volume  (volume  reverberation)  or  near  the  illuminated  surface 
(surface  reverberation).  The  Poisson  assumption  can  clearly  not  account 
for  such  physically  observed  phenomena  as  a  coherent  component  in  the 
specular  direction  ol  a  surface  reflected  wave  or  the  mode  structure 
induced  in  the  scattered  energy  by  a  more  or  less  periodic  wind-driven 
surface.  However,  the  model  appears  not  unreasonable  for  volume 
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reverberation  and  may  adequately  represent  some  of  the  features  of  surface 
reverberation,  at  least  if  the  back-scattered  component  is  of  primary 
interest . 

The  report  deals  at  some  length  with  the  effect  of  signal  bandwidth. 

A  linearly  frequency  modulated  signal  of  wide  bandwidth  (large  frequency 
deviation)  imposes  a  rather  narrow  range  gate  on  the  reverberation.  Its 
width  is  the  correlation  distance  in  the  water  of  the  transmitted  signal. 

Only  scatterers  separated  in  range  from  each  other  or  from  the  target  y  less 
than  this  quantity  contribute  to  the  reverberation  noise  at  the  output  of 
the  cross-correlator.  A  wideband  signal  also  produces  angular  discrimina¬ 
tion  against  false  targets.  This  effect  becomes  significant  when  the 
spacing  between  receivers  is  large  compared  with  the  signal  correlation 
distance . 

When  many  weak  scatterers  but  no  strong  scatterers  (false  targets) 
arc  illuminated,  the  model  leads  to  a  Gaussian  reverberation  process.  The 
detector  output  signal-to-noise  ratio  is  calculated  for  both  reverberation 
limited  and  ambient  noise  limited  environments.  Under  the  stated  assumptions 
(no  Doppler  shift)  the  output  signal-to-noise  ratio  in  a  reverberation 
limited  environment  is  generally  maximized  by  use  of  the  widest  possible 
transmitted  bandwidth.  In  on  ambient  noise  limited  environment,  on  the 
other  hand,  narrow-banu  signals  (no  frequency  modulation)  and  narrow-band 
receivers  produce  the  best  output  signal-to-noise  ratio.  The  output 
signal-to-noise  ratio  varies  as  the  square  of  the  input  signal-to-noise 
for  low  input  signal-to-noise  ratios,  but  as  the  first  power  of  the  input 
signal-to-noise  ratio  for  high  input  signal-to-noise  ratios.  The  latter 
behavior  is  basically  that  of  a  coherent  detector  whereas  ..he  former 
corresponds  more  nearly  to  incoherent  detection.  Even  at  very  high  i: put 
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signal-to-noise  ratios  tin.  correlation  detector  is  still  inferior  to  the 
true  coherent  detector  (replica  correlatoi)  by  3  db  of  output  signal- 
to-noise  ratio.  This  is,  of  course,  due  to  the  availability  in  the 
replica  correlator  of  a  noise  free  replica,  whereas  each  of  the  channels 
of  the  cross-correlator  is  noisy. 

V.  Vertical  Directionality  of  Ambient  Noise 

Report  No.  25  uses  an  ambient  noise  model  introduced  by  Tnlhora 
(JASA  36,  1541,  1964)  to  examine  the  vertical  directionality  of  the  ambient 
noise  in  a  number  of  practically  interesting  situations.  The  model 
attributes  ambient  noise  to  a  series  of  statistically  independent  noise 
sources  uniformly  a' stributed  over  the  ocean  surface.  Talham's  calculated 
noise  distributions  for  hydrophones  located  near  the  bottom  showed  good 
agreement  with  measured  results,  at  least  for  low  sea  states.  Report 
No.  25  is  concerned  primarily  with  ti.e  noise  field  near  the  surface 
(100  -  500  ft.  depth).  Pealistic  velocity  profiles  for  several  locations 
and  seasons  were  considered  The  results  indicate  that  very  substantial 
differences  may  exist  between  the  r>  ' ’-e  intensities  in  the  vertical  and 
horizontal  direction  (differentials  in  excess  of  20  db  were  indicated 
in  a  number  of  cases)  .  The  computed  airectionality  patterns  fall  into  two 
classes.  If  the  observation  point  is  ir.  a  region  of  positive  velocity 
gradient  the  pattern  exhibits  a  peak  near  the  horizontal.  If  the  observa¬ 
tion  point  lies  In  a  region  of  negative  velocity  gradient  there  is  a 
sharp  null  near  the  horizontal.  Except  for  this  anomaly  (which  is  confined 
to  the  immediate  vicinity  of  the  horizontal)  there  are  no  major  qualitative 
differences  between  the  various  calculated  patterns.  Aside  from  the  major 
seasonal  changes  In  velocity  profiles  the  results  do  not  appear  to  be 
critically  dependent  on  the  particular  profile  used. 
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Summary 

This  report  denis  with  a  suboptim-'l  instruru  ntntion  for  the  detection 
of  a  weak  target  in  the  presence  of  interference  iron  a  much  stronger 
target  as  well  as  ambient  noise.  The  instrumentation  considered  is  identical 
to  tie  nulling  detector  proposed  in  iltaort  'Jo.  21  with  the  hydrophone  inputs 
being  modified  by  sampling  and  ;.-'rd-liir.iti  ng. 

It  is  shown  that  this  implementation  improves  detectability  considerably, 

relative  to  the  standard  clipped  power  detector  PCA  ,  for  point  source 

irterference  targets.  The  cost  of  hard-limiting  is  approximately  2  db 

far  an  interference  to  ambient  noise  rower  ratio  ^  of  unity  and  increases 

N 

almost  linearly  with  ^  at  a  rate  of  Z  do  per  It.  cade.  r.nu:>  for  very  strong 
interfering  "j  eint  s  urcos"  ^  JO  db ,  the  cost  cf  clipping  be-acr.es 
significant.  For  this  the  improvement  in  nullirg  relative  to  the  PCA 
is  approximately  29  db  for  a  CO  hydrophone  array.  In  addition  it  is 
demonstrated  that  the  assumption  that  such  a  strong  interfering  target 
can  be  regarded  .as  :  oj  nt  source  is  in  fact  not  unreasonable. 
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I,  Introduction 

Report  K!o.  23  jrrr-c.St.-ri  a  nulling  liti.  ct.^r  Lint.  r- .ximates  Lht. 
likelihood-ratio  d-.tector  for  the  passive  detection  of  a  weak  sanar  target 
in  the  presence  of  ambient  noise  and  interference  from  a  much  stronger 
second  target.  This  rt port  considers  the  detectability  cost  of  sampling 
and  hard-limi  ting  +  he  hydrophone  inputs  prior  to  the  nulling  instrumentation 
Intuition  and  previous  r-. suits  indicate  tint  this  cost  is  essentially 
that  inherent  in  !  -rd-limiti ng.  however,  should  this  cost  be  severe,  it 
will  have  Lo  be  verified  that  the  likelihoou-rntio  test  for  clipped  inputs 
is  nut  sLj.niiicnntiy  better.  Such  verification  is  oxtr.mcly  difficult  ar.d 
will  be  omitted  for  the  present. 

The  nulling  detector  assumes  a  Linear  array  consisting  cf  i-i  equally 
spaced  hydrophon  s.  The  hydrophone  outputs  •  r-  delayed  to  align  the 
interference  components  an-)  or.  then  subtracted  pairwise  to  eliminate  the 
interference.  The  resulting  liffereiicos  are  delayed  one.  nor-,  in  suer,  a 
manner  as  to  align  their  signal  components.  The  outputs  cf  the  s-.  cond  s«.t 
of  delays  are  summed,  squared  and  filtered  in  the  conventional  manner. 

In  the  modified  detector,  the  hydrophone  outputs  ar  sampled  and  nard- 
limitod  prior  to  any  processing. 

Ii.  Analog  hulling  3  t.—tor 

The  detector  proposed  ir.  Report  Wo.  21,  which  will  be  called  the 
Analog  lulling  Ditoctnr  (D  ),  calculates  the  fcl  lowing  test  statistic 

:rs  |  n-r  '  2 

)  3  |  _  x.(tHko\  , 
k=l  I  i=3  , 

where  T  is  the  decision  time,  f  is  the  sampling  r-'te,  o  is  thk.  sample 
interval  (  -  -  1  / f _  ) ,  and  i-i  is  the  number  of  hydrophones  in  the  array. 
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The  processes  x^(t)  -ire  given  by 

xi(t+k0)  =  £  vxco  -  \  yijr(k)  , 

where  y^(k  )  represents  the  out  'ut  of  the  i^1  hydrophone  ielayed  first  by 

ia,  where  a  represents  the  delay  per  hydrophone  needed  to  .align  the  array 

with  the  interference,  and  then  delayed  again  by  iX,  where  X  represents 

the  delay  ncr  hydrophone  needed  to  align  the  new  array  with  the  signal. 

The  process  y.  (k)  represents  the  output  of  the  (i+r)^  hydrophone 
^  > r 

delayed  first  by  (itr)c  and  then  •’gain  by  iX.  Hence 

yi(k  )  =  i  ^  ( t  v  kc  +  iX  +  iu )  +  i  ( t  +  k  a  +■  i  X )  *•  s(t  +  ko)  , 

nnd 


y.  (k)  =  n.  |t  ko  +  iX  *•  (i  •  r)a’«-  i(t  t-  ko  +  iX)  +  s(t  +  ko  -  rX; 
i,r  i  *-  r  j 

where  n(t )  represents  the  ambient  noise,  i(t)  repr.  s  nts  the  interference, 
and  s(t )  represents  the  signal.  If  either  the  interference  or  signal  ware 
broadside,  tne  second  delay  would  be 

X  =  —  sin  -J  , 
c  ’ 


where  d  is  the  sr icings  between  hydrophones,  c  is  the  sored  of  sound  in 
the  ocean,  -'r.d  Q  is  thi  angle  between  the  target  and  interference.  The 
procosse  n  (t)  v.  ,r  -  1 ,  i(t);  and  s(t)  are  assumed  to  be 

statistically  independent,  normal  processes  with  normalized  correlation 
functions  p  (  ),  p^(  ),  and  p  f  )  and  whose  po  .-rs  r>re:  N,  I,  and  3  respectively 
Let  us  first  develone  an  expression  for  the  autocorrelation  function 


x.(t*kcj  defined  as  ti.  .(k). 


It  is  seen  that 
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Ri  =  i  '■  yiTo^FXk)  -  y,  (o)y  ( k )  -  y.  (o)y  (k)  *  y.  (o)y  (k)  ,  (5) 

4  i  1  •'  1  1>4  J  4Jr  J)v  i 


where 


y~ToJy  (k)  =  N&(i*n-j-n)  ,■»  ko+(  j-i  )X+  (  j+n-i  -m  )a'+lp.  f  kc+(  j-i  h+Sp  |  k*-+  (ra-n  )\i, 
>  ITl  J  .?  11  11  lit  -  Sl 

we  can  distinguish  three  cases  of  interest. 


Case  1  (l  -  j ) 

Ky_(k)  =  \  '  ilp„( ko )  +  S|p<3(ke)  -  \  (ko+rXj-  \  Pc(kc-rX) 


2  ''n 

i. 


2  -  s' 


2  Ks' 


r  • 


Case  2  (  i  +  r  =  j ) 


»<i  f(k)  -  |  3 '  p  ^  ( ko )  -  |  ps(  knt-rX )  -  ~  p.(kc-rX)1  -  |  .:pn(kc+  i\)  '  , 

Case  3  ( ,1  -  i  +.  p;  p  n  r  ) 

1  at 

,  .  D  (k)  =  ^  5  p  ( kn)  -  ~  p  (ko+rX)  -  x  p  (ko-rX)  . 

2i  £  b  2s  i 

It  is  convenient  to  -jejune  R.  .  (k)  =  it.  .  (k)  -  R  .  (k),  where 

1,1+ r  -  1,1+ r  a  ,i+_  p 


(7) 


(d) 


(9) 


(10 


R.  .  (k)  -  -  r  N  p  (ko  +  rX) 

l ,  i  *_r  h  n  - 

ue  nr-  now  in  a  position  to  evaluate  the  si pnal-to-noise  ratio 

SNR  af  the  r.ullinr;  o*  t-ctor  for  weak  input  signals.  This  ratio  is  defined  by 

-  B  *  '-11} 


where 


(.  =  1  im  3fx  )  -  E  ' 3 \  x  )J  , 

|  ,  ,  r.  ,  , 


'ind 


'jruJ0'(w  . 


"lie  subscripts  K  and  h  r  for  to  the  hy-otN  sis  of  ,i,  nal  ..b.:ont  t  i.d  the 
•alternative  of  signal  rci  It  is  easily  s*-on  that 


M 


x  tTx  .  (  t  )  '  -r 


Tf 


C.,  Z_ 


M  -r 

r 


(ii) 


UO 


)R  .  ( c  )+2  (1-i-r-p  )R.  .  (o)+2(H-2r)R  .  (o). 

a,1  4^.  a?t+p  i)1!1 


(lb) 
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Hence  substituting  the  r.  salts  ci'  iou  -t-or.s  •«< J  10  into  p:>t  Lon  lit 

and  then  inti  Tquit Lor;  12  /iells 


(  r  )2  ' 

A  -  If  -v— 5^  1  -  p  ( r>  ) 

an  s  2  ,  .s  I 


(19  ) 


The  variance  of  T  (x)  is  riven  by 
nn'  ‘  J 

Tf 


-.r  is  (x)l  -V  V 

I  an  i_ 

k  q 


,  p(k  -  |Hk)2JHq/ 


(16) 


where 


P  (  k  ) 


X1 i t t Kg  ) 


i=l 


5ir.ee  fl (k)  is  i  nn  i,  L  •  follow.*  ti.st, 

-  r 


■  V  1  6  [  x  / 1 


k  q 


r  ( K  )|Kq )‘ 


Tf 

_o  s 


(k) 


Pf 


P=1 


t’(o)i3(r  )  •  ,  (i-  pf-'  '> 

p=i 


(17  ) 


Wh  .IV 


p(0  ;p( S  ■ 


.  .-r 


i  J 


i  .-r 


x.  { t  x  .( ••■-k  •  -  ( ! '  —  r  ill.  .  (si 

i  t  ,  . 


(k-r-r  .  {, 

i  ,  i-  p 


o  =  l 


2;.-.-2r,h.  . 

i  -r 


!  ]  M  5 


Once  torn  cc  oi  r  ir.r  ..  i.uatior.s  '.',>,10,13,17,  ar.'i  Id  results  ir., 

Tf 

2  „  f,  _  .2  s  k  ; 


.2 


..-r  .  ,2  .  ;  ,  ,  i i — 2 r 


: i  ■  r  4  ’  /  ‘  •  *“  i  > 

•  rl3(  — , 


I  L 


1-  (  flF  )0rS  rX }  -  2  •  (  "Tf  •  pr.(  kT  } 


k=l 


(1?) 
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From  Equations  11,  l1?,  and  19  on..  fin  illy  obt;  irs, 


Tf  „ 


(  -r)fl-p  (rX) 


Tf„ 


(rX)!  -2 

^  i-.  -  r  n 


(1- 


Tf 


•)'P  (kn)-  i([— — )Jp  (ko+rX)+p  (ko-rX)f 
'  'n  2  Ii-r  |  n  n  ( 


1 1 


l'  k=l  •  L  1  (2to) 

Throughout  this  report  we  will  m,ake  the  .assumption  that  rX  is  sufficiently 
largo  so  that  Pn(rX)  p^(rX)  =  0.  This  is  not  the  same  as  requiring  that 
the  target  and  interference  ar<  f  t  r  ncvt  i  in  bearing.  ,-cr  tt.is  later 
condition  ths  assumption  is  satisfied  for  r  =  1,  and  wc  can  subtract  a  g  '-cent 
hydrophones.  When  tt.e  interference  is  relatively  close  to  the  target,  the 
•assumption  can  be  satisfied  by  making  r  large  enough.  For  large  arrays 
this  can  be  accomplished  at  a  fairly  small  cost  for  all  bearing  angles 


except  those  in  the  immediate  vicinity  of  si.ro.  For  example  if  the  it.t  uts 

are  pre -whitened  arid  processed  ut  to  >  KC,  tlic  hydrophone  spac-ings  are 

2  feet,  and  r  =  h;  the  assumption  is  identically  satisfied  for  a  relative 
o  „  I 

bearing  angle  of  3  3v  in  the  vicinity  af  broadside. 

The  primary  motivation  for  the  previous  as,' a unr ti on  is  that  it  will 
greatly  simplify  the  analysis,  particularly  when  we  considi  r  hard-limited 
data.  The  justificat  on  can  be  found  in  the  discussion  of  Report  no.  21 
■about  the  bearing  response  pattern.  It  was  pointed  out  that  while  the 
SNR  of  Equation  20  is  insensitive  to  the  relative  bearing,  the  baring 
response  pattern  Iooks  rather  had  for  r  =  1  when  the  interference  is  relatively 
close  to  the  target.  It  was  further  noted  that  the  hearing  rcspcr.se  can 
be  greatly  improved  by  increasing  r.  In  practice  the  bearing  response  is 
suliiciently  important  as  to  require  the  assumption  that  pn(rX)  '=  Ps( rX)  =  0 
for  all  bearing  angles  above  seme  small  value.  For  this  implementation, 
Equation  2d  car.  be-  closely  approximated  by. 
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SNR  '*  ~  % 
an  I  2  M 


I  T  f  ,,  '1 

!,  ,  'S/,  k  2,  ,  l,i-.-2r  2,,  -.2 

i  t  2_  (1-  — ),Pn  (ko)+  2  ,TC~r  pn  (ko~rX)  . 

I  k=l  3  1  1 


Let  us  new  consider  the  special  case  where  tne  inputs  are  pre-whitened 
to  wq  radi:  nr/sec  •  r.d  sampled  at  the  Nyquist  rate  (f  =  w^/r. ).  Equation 
2}  becomes 


,Tw 

“  =  J—2.  ± 

an  r  2n  N  r 


_ l-i  -  r _ 

2  1 
,  1  ,h-2r  v  /,  r\  .  77 

1  +  2  (IT07)  (1  "  —  )  - 


and  for  targe  decision  bira  s 


—  ^  \  5  M  -  : 

Sink  «  '•  — —  —  - rr 

3n  1  cTT  IM  /  -1 


/  1  ,M-2r,  2' 

A  f  2  (fTT7' 


Replacing  r  by  1  results  in  Equation  20  of  Report  No.  21. 

As  a  second  case,  lot  us  rewrite  Equation  21  for  large  decision  tir.es 
in  the  f  11 owing  way 


L>.  ft  = 
an 


W  j  (K  -  r) 


‘l2o  +  V_  [t.‘<l;")0  + ! 

L  k=l 


M-2r|2  2 

l~)  P" 


(ko-r A)  o 


For  extremely  rapid  3 am pi Lng^rates  ( <--s  .entially  continuous  operation) 
this  approaches 

t  M  _  r 


an  1  or 


a  on  (  p)dp 


-  /,  ,M.-2r  N 

V1  4  (*r?> 
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It  will  become  a... parent  th  >t  even  for  pre -whitened  inputs  one  should  sample 
significantly  faster  than  the  lyquist  rate  for  large  interference  to  ambient 
noise  power  ratios. 

III.  Clipped  hulling  D  tector 

The  modified  nulling  detector  will  be  called  the  Clipped  !  ulling 
Detector  (D  ).  The  first  modification  of  the  previous  analysis  occurs 
at  '  u--ti.cn  6.  Since  the  proc  sses  y(k )  are  limitod  to  the  values  1  or  0 
depending,  on  tfu  ir  sign,  Equation  6  becomes 


~  (n  )y  .  (  k  ) 

x ,  m  .1 ,  n v 


2  -^r N  ( i+m-j-r.  jp  ko<-(  ,i-i  )X+v .iHi-i-n  )aKp.  .<e+  j-i  )X+Sp= ko+(  rn— n  ) 

-si n”x  J - ™ - - - — - — - - - - 

L  N  +  I  +  S  (2I4) 


It  follows  that  for  the  thr<_e  cases  of  int  1  st 


ase  1  (i  =  ,j) 


,  sin 


r,.-i 


x  ,  ! 


•  !p  (ko  )+ip.  (  kc  )+Sp  1  ko  )1 


.-i  +  I  +  3 


|  sin'1 


Ipi(ko)+Spa(ko-r\) 
.1  +  I  +  S 


1  .  -1 

-  2  sin 


Ip .  (  Ko  )+  Sp  (  kc+r  X)a 


iJ  +  I  +  o 


Mf  • 


(25) 


Case  2(,j  *  i  +  p ;  p^r) 


R.  .  (k)  =  - 

1,1+ p  n 


_a 


1  s  in 


Ioi(kc  +  pX)+Sps(  kc  )|  x  _1j-Ipi(ka+pX)+Sos(kc-;TrX)-j 


N  -  I  +  S 


-  2  sin 


N  +  I  +  S 


1  .  -1 
-  2  Sln 


Ipi(kc+pX)+3f>s(kc+rX)-j  | 


ase  3 


(,i  =  i  +  r)j  R 

L}1 


i  +  r 


H  +■  1  +  O 

(k)  =  R.  .  ( k  )  -  R .  .  (k) 


IJ  . 


(26) 


1,1+r'  1,1 +p 


1 

2m" 


r.  -1 

{sin 

rlp^(ko+rX)+3ps(ko+r\)‘ 

.  -1 
_  ei  n 

L  N  +  I+o 

“  Jill 

L 

ilpn(  kc+ rX  )+IPi(  ko+rXj+S^knr  >'■]  I 


N  +  I  +  S 


(27) 


Subst j tut-*  ng  these  three  Equations  into  Equations  12  and  ].U  yields 
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A  =  Tf  — 
cn  s  n 


I  .  -1 

sin 

L 


N+I 


.  -1 
-sin 


I+5P3(rX)il  ,M-2r J  __._-irI^n(r>-)+Ipi(rX)l 
N+I+S  Jj+  M  -  r  j Sin  |_  N  +  I  J 


-  sin 


sin 


_lfN|»n(r'O+lpi(r\)+Spg(r>0 


iPi(r\) 


il  +  1  +  J 
i-i-r 


.  -1 
+  sm 


flp  (rX)+Sp  (rX)a 

1  O  . 


N  +  I  +  S 


,  + 1  ir2I.(1 

J  k=l 


Ipi(kX)+S 


N+I+S 


.  -1 
-  sin 


Ipi(kX)+Sp3(rX) 


IJ  +  I  +  S 


(2^) 


Making  the  assumption  that  p  (r\)  =  p  (rX)  =  p.^r\)-0,  simplifies  the 

ns?. 

previous  equation  to 


i-i-r 

r 

,'sin'1 

[  I 

.  -i 

,  M-r 

1  H+pVti  k 

.  -1 

Ip4  (k\  )  +  3i 

Ipi('-X)  - 

n 

N+I 

-sin 

N+I+S  '  r 

sm 

A  -  I  «-sl  Sln 

>  j  +  s 

/ 

J" 

* 

L  '  k-1 

L  .J 

(29 

y 

if 


It  is  readily  shown,  via  the  Taylor  expansion,  that 

A 


.  -1  A+bs  .  -1  A  \ 

Sln  N  +  I+S  =  Sln  %+l) 


b  - 


N+I 


=?+  0vN+I^ 


Hence  Equation  29  can  be  further  simplified  to 

ti-r 


a  t -r  I  *  *-  r  S 

Acn  "  s  n  N+Il 


J/i^1  +  2V  (i  _  JL.)  -■  ■■ — 

"Vl— (I/N+I  )2'  ^  h_r  1  7 


k-x  jj^fikx)'. 


(30) 


This  can  be  more  conveniently  approximated  by 


cn 


=  Tf  —  ~  J  rt-r-l  +  I/N  +  I 

s  n  N+I 


K-r 

I 


cx  .  k  pi  (kx) 

,/  fT/M  T  n2'  ^fl+I  ^  2_,  ^  ~  FUr  ^  r  t  "  - - 1  I 

yi-(  I/N+I )  k»l  j)2Pi2(kA  T  J 


(31) 


We  now  wish  to  calculate  the  variance  of  the  test  statistic  Scn(x). 


Unfortunately  the  nonlinear  properties  of  hard-limiting  make  this  calculation 
for  the  general  case  impossible  to  perform  exactly.  However  by  suitable 
approximations  we  can  obtain  a  reasonable  estimate  of  this  variance.  Let 
us  first  consider  a  very  special  case. 
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Pre-Whitcned  Inputs  ari  Slow  -sampling  Rites 

The  assumption  of  p re-whitened  inputs  out  to  the  processing  frequency 
is  certainly  reasonable  if  one  is  concerned  with  the  inherent  cost  of  a 
particul'  r  implementation.  Let  us  furtner  assume  that  the  inputs  arc- 
sampled  at  the  Nyquist  rate,  chat  is  all  of  the  samples  in  the  output  of 
a  single  hydrophone  ar<  statistically  independent.  Since  the  sampling  r  -to 
determines  the  minimum  delay  it  quantizes  the  permisf-ablc  bearing  angles 
i.o  >•:;  1  oh  tie  rray  c'u  be  aligned.  Thus  point  sources  can  be  only  cru  K.'iy 
rulli  d  oet  even  in  the  absence  of  hard-limiting.  It  will  be  se_n  that  for 
interf-  ren.es  to  ambient  noise  ratios  significantly  greater  than  unity, 
this  "sampling  cost"  severely  degrades  the  performance  of  the  det-ctors.  It 
is  useful,  however,  to  consider  this  restrictive  case  first.  Fina  lly  let 
us  -assume  that  the  array  can  be  aimed  directly  at  the  int-.r rcr^nce.  This 
a  isumption  not  only  removes  the  "sampling  cost"  but  greatly  simplifies  the- 
analysis . 


Under  the  hypothesis  (3  =0),  given  the  assumptions,  equations  23, 
26,  and  27  can  b .  rewritten  -s 


Ri,i(k)=  I 
Ri,i+p<k>  =  °> 


R.  .  (k)=  <t 

i,i+r  1 


a  2  .  -1  I 

1  -  -  Sin  -r—r 
n  r!  t- 1 


1 

l 


,  2  .  -1  I 

1  -  -  sin  — p 
n  H  +  I 


d  . 


where  X  =  mo.  Since  \  S-  -  it  follows  that  m  ^ 


k  =  0 
otherwise , 

if  k  =  +  rm 

otherwise, 

d 


(32) 


co 


The  variance  of  3 


cn 


under  the  hypothesis  is 
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A1#-  ,(l,-r,) 

VarH^Scn^=Z!Z^  Y*)  ),  }  Xj  (t+ko)Xj(t-*-ko)xr(t-t-qc)xs(t+qg)-x1(t  )xj(t )  xr(t)xg(t)  . 


k  q  i  j  r  s 

Consider  first  the  term  t^,  obtained  by  setting  k  =  q. 
(M-r ) 


(33) 


x  (t  )x  (t  )x  (t  )x  ( t )  -  x  (t)x  (t)  x  (t)x  (t) 
ijrs  ij  rsj 


From  the  feet  that  x^(t)  is  either  1,  0,  or  -1  and  from  Equation  32  we 
see  that  x^2(t)  =  1  with  probability  p_=R  (o)  =  ^|l-  —  sin  1  yyy  ]  and 
atro  oti."  rwise.  It  follows  that  (t)  =  1  with  probability  p.  Hence, 


yy  =  (:\-r  )(p-p2  )  +  2^^xi^(t)  xs2(t) 
=  f2(n-r)2-  3(H-r)lc2+  (M-r)p. 


(3b) 


l’he  term  t^  obtained  by  seating  k  ^  q  and  ignoring  the  special  case  where 
k  -  q  =  j_  rm  can  easily  be  shown  to  be  identically  zero.  For  this  special 

case,  1  t.  us  consider  q  =  k  +  rm  and  by  symmetry 

_ (M-r) 

v  ( t  )x  t  )x(t+rmo  )xg(  t  +  rmo )  -  xi(t)x_.(t)  xr(t)xg(t)j  . 

i  j  r  s 

This  term  has  a  nonzero  value  for  r  =  i  -  r,  s  =  j  -  r,  and  i  /  j  (s  and  r 

can  be  reversed)  as  well  as  for  the  case  i  =  j  and  r  -  s  =  i  -  r.  Hence 


t  =  U(  Tf  -rm ) )  )  x.(tjx.  (t+rmo)  x.(t)x  ^(t+rmo' 

j  s  L — <Z <  i  i  J  J"’1'* 

i/ J 


2(ffs-rm)  y>i2(t)x._r2(t.rmo)  -  x±2(  t  )xi_i,2(  t ) 


recognizing  that  R.  .  (rm)  =  -  r  p,  we  obtain 
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t3  -  (Tfs-rm)[(h-Lr)2-  (h-2r)]p2+  (2Tfg-rm  )(li-2r  ^  , 


(35) 


where 


'L  =  \  (t)x._r  (t+rma)  -  (t)  xi_r"(t)  . 

Combining  Equations  and  35  we  obtain 


Dc„2-  2T‘A  -”)2p(p 


1- 


rm  wII-2r 


21 


2  (M-r)  2U“  Tf  ^  M-  r  } 


+  ^N-r) 


(36)' 


Although  we  cannot  eval  itate  we  do  recognize  that  2p.  -p11"  is  quite 

o 

small  (of  the  order  i  ^  T  In  fact  if  x(t)  were  normal  (not  true)  the  torn 

2  2  1 
2ll  _p  would  be  identically  zero.  Finally  since  ^ ^  is  noriv^lly  much  less 

I  2 

than  p  (except  for  very  large  —  in  which  case  2M—  p  is  bound  to  oe  very  much 

2 

less  than  p),  the  error  ir  assuming  that  2m— p  =  0  is  quite  negligible.  There¬ 
fore  to  a  very  good  approximation 


D 


cn 


|(1-  yXfrn:!2*  • 


(37) 


Combining  Equations  11,  31,  and  37  we  obtain  that  for  the  ;  re-wnitened  case 
with  slew  sampling  rates 


SNR 


2  R  S 


M-r-l+ 


l A',*  I 


cn  n  V  2n  M 


(1*J)0-  -Sin"1  rri-) 
N  n  b  +  I 


,  l,M-2r,'  ,,  rXs  l-3p 

1+  2^T3T)  ^1-  t-) 


(38) 


?(M  -r)p 


1  2  IT 

rfhere  p  =  ^(l-  ~sin  'r7x')*  The  cost  of  clipping  for  these  slow  sampling  rates 
Is  obtained  by  comparing  equations  38  and  22.  If  the  cost  (C  )  is  defined 
as  ti  e  decrease  in  SFR, 

I/N  +  I 

h-r+  j— 

,  2  1 


/ "  Iw,  2  .  -1  IT 

(1VHl"  n  3in  ii+T 


/ - -  /.  l/M-2r. 

Vl-(I/N  +  r  )-.  /  1+  2  i,  -  r  j 


M-r 


4+l(li=i£)2+_l=3E_ 

2  rt  -  r  2  2 (M-r  )p 


(39) 
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For  large  arrays,  equation  i9  can  be  approximated  by 

1 


l 


"  <1+  ffX1  "  \  sin_1  HTT> 


(Uo) 


If  we  now  assume  that  —  is  large  and  consider  only  the  first  two  tor  s 
of  the  expansion 


co 


.  -1/-  \  n  T-1  f,  1-3*5  (2k-l)  _k 

sin  (x-a)  =  -k  -V2z  1  +  >  -57- - ^ - L  z 

d  I  — *  2liic(2k+l)  k.i 


k=l 


we  obtain 

Cm  —i -  .  (la) 

V^TTl) 

It  is  observed  that  the  cost  of  clipping  increases  with  i  and  gets  quite  large 

for  large  values  of  /bile  this  cost  will  get  cut  roughly  in  half  with  more 

sensible  sampling  rates,  it  will  remain  significant.  The  following  table  is 

an  evaluation  of  Equation  39  for  a  variety  of  values  of  ^  and  M  under  the 

assumption  that  r  =  1.  (Actually''  the  assumptions  inherent  in  Equation  39 

make  this  choice  the  only  logical  one,  since  the  sampling  rate  prohibits  the 

detector  from  searching  for  a  target  at  any  bearing  angle  other  than  those  for 

which  p  (X)  =  0.  Eouations  38  and  39  are  of  course  based  on  the  assumption 
n  — 

that  any  target  present  is  precisely  at  one  of  these  bearing  angles.  These 
restrictions  will  soon  be  relaxed.) 


0 

1 

10 

100 

100C 

10 

1.9 

2.9 

6.15 

10.2 

13.9 

Uo 

2.0 

3.15 

6.6 

11.1 

15.5 

100 

2.0 

3.2 

6. 7 

11.  U 

16.1 

2.0 

3.2 

6„? 

11.5 

16.5 

Table  1  Cost  of  Clipping  in  db  for  Slow  Sampling  Rates 
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While  we  are  considering  slow  sampling  rates,  let  us  calculate  the 
the  improvement  attained  by  nullini,  relative  to  the  standard  clipped  power 
detector,  the  polarity  coincidence  array  (D  ), 

P  C  3 

IV.  Improvement  due  to  Nulling  at  Slow  oampling  Rates 
The  PCA  calculates  the  statistic 


Tfsr  M 
k=l L  i=l 


(U2) 


where 


and 


x±(t)  =  sgn|yi(t  )|  , 


y^(t)  =  n^(t+ia)  +  i(t+iX)  +  s(t)  , 
and  where  a]l  the  symbols  are  as  before.  Hence  it  follows  that 


x.  (t  )x  .(t+ko  )  =>  -  sin-1 
i  J  n 

It  is  readily  Seen  that 


fW6i/n(ko)  .  IP, 

ko+( j-i ) X ] 

+  Spg(ko)  ' 

[  .1  +  I  +  S 

E(3  )  "  _  " 

Trf2-  ’  j(l)  "  H  * 

i  j  k=l 


Jl.  _,rlp  (kX)  +  S 

nL(rf"k)  Sin  '  nVT";~S 


Hence  from  Equation  12 

M 

A  =  lim  Tf  —  (h-k)Jsin  ^ 
Pca  3/N— )  0  5  > 


Tp^(kX)  +  S 


iT  +  I  +  S 


sin 


-1 


flp  (kX) 


N  +  I 


(a  3) 


(Ul*) 


(W) 


It  follows  that  for  pre-whitened  inputs  and  f'yquist  rate  sampling  that 

A 


=  Tf  ~  ~  (M2-  M) 
pca  s  7i  N+I 


(U6) 


The  variance  of  S^ca  under  the  hypothesis  (S  =  ij  )  and  v;ith  the  slow 
sampling  assumption  is 
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Var„ Is _ ,(x)"|  =  )  )  )  )  )  /  x.  (t+ko)x  ,(t+ko)x  (t+qo)x  (t+qo);-  M2 

n  i.  j  4 — ,  t ,  £ i  l — i  / —  c — i  i  j  r  s 

L  i  j  r  s 


where 


and 


i  ;  i  t  j 

X  (t  )x  (t  ) 

3  1  o  ;  i  /  j 


x^t  )x  (t+j>o) 


2  .  -1  I 

-  Sin  rr~7 
n  N+I 


j  =  i  +  p/m 
otherwise 


Since  x.  ’t)  =  1  for  all  i,  equation  li?  can  be  rewritten  as 


.  ,  ,  ,  V  7  x.(t+ko)x  ( t+qo  )x  .( t+ko  )x  (t+qo) 

p°a  /L/L  C,L>l,L  1  r  y  c 

q  l/j  rfa 


2Tfg^7x7  (t)  xT  (t)  -*■  l^(Tfg-pm)^)  £x.(t)x._p(t+pma)x.(t 

p=i 


2Tf 


2  q 

(M2-  M)  1  +  2(  --  sin”1—)  £(1-  *p)(M-p)(M-p-l) 

P-1  J 


Recognizing  that  2 y  ( M-p  ) ( ii-p-1 )  =  ^ fffilki  (H^-  m),  Equation  Li9  becomes 
P-1 


large  decision  times 


D  =  2Tf 
pea  s 


1  +  M  <  !  sin_1  nTI)2 


2 

M  -M 

where 

M  =  (2M  -  U)/3  • 

finally  from  Equations  U6  and  50  we  obtain 

V  M2-  M _ 


SNR 


2-\/Tw0  S 


pea  7i  V  2n  N  +  I 


->/7 t  /2  7-1  I  \2 

V  1  +  -1  n  Sln  N+ 
A-L-l 


,  (U7) 


(U8) 


)x^(fr-pmo) 

(U  9) 


for 


(50) 


(51) 
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Comparing  this  res  lit  with  Equation  j8  we  see  that  the  improvement  as  a 
consequence  of  nulling  (j)  for  slow  sampling  rates  is  given  by 


i 

Q  - 

i- 


77  ,2 


1  +  M  (—  sin 
n 


-1 


I/N+I ) 


I/N+I 


~v£(i7tt+i) 


T\ 


1  +  1/M=2 £)2+  l-3p_ 
T-  M-  r '  2 ( M-r  )p 


‘ \ 


V  M2- 


M 


,  2  .  -1  I 

1  -  -  sin 

n  N+I 


This  improvement  is  plotted  in  db  in  Table  2  for  r  =  1. 


0 

1 

10 

100 

1000 

10 

-1.6 

1.6 

7.8 

lli.2 

20.8 

hO 

-1.0 

3.7 

10.6 

16.6 

22.  U 

"00 

L. 

-.8 

S.h 

12.6 

13.5 

23.9 

(52) 


Table  2  Nulliru  Improvement  in  db  for  Slow  Sampling  Rates 
V.  Effect  of  Sampling  Rata  on  Nulling  Procedure 

As  mentioned  before  the  sampling  rate  quantizes  the  bearing  angle  to 
which  the  array  can  be  aligned.  Let  us  therefore  assume  that  the  array  was 
not  aligned  exactly  with  the  interfering  point  source.  In  the  analog  nulling 
detector  the  difference  of  a  pair  of  hydrophones  is  seen  to  have  the  power 

x2(t)  =  2n|i  +  -  Pj  (re  )]|  >  (53) 

where  e  represents  the  error  in  the  delay  per  hydrophone.  If  a  is  the  sampling 
interval,  it  is  clear  that  the  maximum  error  in  the  delay  between  adjacent 
hydrophones  is  o/2.  oince,  in  principle  the  hydrophones  need  not  be  all 
aligned  in  one  direction,  the  maximum  error  in  the  delay  for  hydrophones 
further  apart  is  not  obvi "US.  We  will  bypass  this  difficulty  by  assuming 
that  the  interference  and  target  are  sufficiently  removed  in  bearing  ..hat  r 
can  be  set  equal  to  one.  Now  assuming  that  the  interference  spectrum  is 
flat  out  to  the  processing  frequency,  we  have  that 
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x  (t)  <:  2N^1  +  i 


..  sin  w0o/2 

1  -  Uqo/2  ' 


]  . 
J 


(5U) 


Previously  we  have  assumed  sampling  at  the  llyquist  rate  (o  =  -^-).  Hence 

w 

o 

the  power  of  x( t )  can  be  as  much  as  (l  +  .  363X  I/N )  times  greater  than  the 
idealized  case.  Therefore,  for  interference  to  ambient  noise  ratios  much 
greater  than  unity,  the  interference  will  still  be  dominant. 

A  more  reasonable  sampling  rate  is  one  for  which  the  "sampling  cost" 
is  at  most  3  db  (when  c  =  o/2).  Actuaxiy  this  sampling  cost  is  not  only  a 
function  of  the  increased  power,  but  the  residual  interference  introduces  a 
nonzero  autocorrelation  between  the  hydrophones  as  seen  for  example  in  equation 
h3.  This  autocorrelation  "costs"  an  additional  amount  by  increasing  the 
test  statistic  variance  (see  equation  30*  for  example).  In  fact  an  increase 
in  the  power  of  x( t )  of  50  per  cent  will  "cost"  in  the  order  of  3  db.  For 
the  sake  of  simplicity  we  will  assume  that  a  reasonable  sampling  rate  is  one 
for  which  the  power  in  the  modified  channels  x(t)  is  at  most  50  per  cent 
larger  than  the  idealized  case.  It  follows  that  a  minimum  reasonable  sampling 


rate  f  is  riven  by 
sm  J 


1  - 


sin  w  /2f  n 

o  sm  1  N 

w  /2f  =  2  I 


(55) 


For  large  I/N  ratios  thus  equation  can  be  approximated  by 

,2 


or 


w 


2f 

sm 

?nf 


sm 


w 


AJ  N 
J  I  • 


~  -,/TT 

v  3  I 


(56) 


These  results  are  tabulated  below 


I 

N 

1 

10 

100 

1000 

2nf. 

oFH 

1.7 

5.7 

18 

57 

Fable  3  "Minimum"  Sampling  Kates  for  Nulling  Procedure 
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It  follows  that  for  an  intf rfc rence  Dower  20  db  greater  than  ambient  noise, 
a  sampling  rate  of  the  order  nine  times  che  Lyquist  rate  is  called  for. 

Let  us  now  consider  the  sat  cling  cost  when  hard-limiting  is  introduced. 


The  equivalent  to  liquation  CU  for  D  is 

cn 


x  (t)  <  2<  1 


1 


sin  w  g/2 
o 


2-,  —  J  lj 

„ .  -1  I  _ o 

n  Sln  N+I  w  o/2 
L  o 


"his  can  be  closely  approximated  by 

T77 ;  _  X  2  .  -i  i  2  Sln  wo°//2 1  i./  ii  ' 
x  (t)  <  2<1  -  —  s:n  —  +  -  1  -  .  ^  — 

If  we  sample  at  the  rate  determined  by  liquation  tnis  becomes 

2,  .  .  _L  2  .  -1  I  1  fN  1  , 

x  ( t  )  <  2<  1 - sin  r;— T  +  ~  (• 

I  n  N  +  I  n  *  N+2I 

Using  the  expansion  lor  5in~\l-e  ),we  see  that  for  large  I/N, 


^  olL  J  N  '  1  /~N  1 

1  ^  '  On  V2N+  21  +  n*V  h+2I  ’  * 


Thus  the  maximum  "sampling  cost"  of  the  clipped  nulling  detector  when  sampling 

at  the  rate  f  is  approximately  1/2  that  of  the  analog  nulling  detector. 

The  implication  of  this  result  will  be  clear  later. 

VI.  Cost  oi  :.ard-Limitinp  with  kapid  sampling  Hates 

The  oroblem  is  to  evaluate  the  variance  of  the  test  statistic  S  (x) 

cn 


under  the  hypothesis  S 


-•s  r’entioned  Before,  an  approximation  is  necessary 


to  overcome  the  analytic  difficulties  of  hard-limiting.  Tne  test  statistic 
Scn(x)  calculates  the  sum  of  the  squares  of  0(k)  where 

M-r 

R(k  )  -  ^  xi(t+ko)  ,  ( 

i=l 

and  where  y.^(t)  is  the  difference  of  a  pair  of  clipped  hydrophone  outputs. 
Under  the  hypothesis  S  =  0,  the  autocorrelation  function  of  x^(t)  with 
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x^(t+Ko)  is  determined  by 


R.  . ( k  )  =  —(sin 
1,1  n ' 


-1  Npn^kc^  +  IP^ko)  _1 

-  -  sin" 


1;  +  I 


flpi(kc ) 
N  +  I 


(63) 


,i+p 


(k)  =  0  , 


R.  ,,  J  k)  = 


if*  _!  [■  kPp(  ko^rX)  v  Ipi(ka+rX)'j 


2n 


i 


sin 


.  -1 
-am 


Ip^ Ko+rX) 


M  +  I 


N  +  I 

where  R.  .  (k)  =  R.  .  (k)  -  R.  .  (k).  It  is  clear  that  the  process  x.(t) 

l , l^r  i,  l  +  r  1 ,  i +_p  '  .  i 

has  an  a.  pliti.de  density  function  given  by  Figure  1,  where  the  amplitude  of 

4\ 


1 


Figure  1  amplitude  Density  cf  x^(t) 

the  side  delta  functions  is  1/2  n.  .(0).  It  seems  reasonable  to  assume  therefore 
that  for  large  arrays  (M-r  >->  1),  p(k)  is  approximately  normally  distributed 
with  zero  mean,  '-ie  will  see  that  for  slow  sampling  rates ,  tne  results 
obtained  with  this  assumption  are  similar  to  that  of  equation  3?  except 
that  the  last  term  is  missing.  ihile  this  last  term  cf  Equation  3?  does 
indeed  vanish  for  large  arrays,  it  can  be  seer.,  by  putting  in  nu.ibers,  that 
for  very  large  T/N  ratios  (30  db )  and  reasonable  sized  arrays  (to  =  hO), 
that  it  is  not  negligible.  This  can  be  traced  to  the  fact  that  the  amplitude 
of  the  side  delta  factions  of  Figure  1  gets  extruely  small  for  large  I/N 
ratios.  ’7e  will  ignore  this  problem  for  the  i  oment  since  it  will  soon  turn 
out  that  this  difficulty  can  be  patened  up. 

The  autocorrelation  of  0(1^)  with  ^(i^+k)  is  determined  with  the  help 
of  Equation  60  to  be 

i-i-r 

p(o)  p(+k )  =  ^  ^  xi(t  )xj(ttko^ 

i  3 
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,j_r  [  -if  N<VKo)  "  lfVk°) 
“—.sin 

’  L 

ii-2r 
2n 


^  sir 
L 


.1  «•  I 


-i  I  Ipi^: '"rXn 


.  -1 
sm 


—  pi(^} 


r  I”2 - 

l  K  *  I 


-  sm 


-1 


Npn(kc-rX)  +  Ip^(ka-rX) 


N  +  I 


Since  p^  (rX)  is  assumed  to  be  negligible. 


P2(k)  «  ‘-^1  | 


2.-11 
-  -  sin  . — - 

Tl  W  +  I 


(61) 


(62) 


we  are  attempting  *  >  u.:  rr.in.  variance  cf  3  (x)  under  the  hypot’i  sis>o; 

cn 


t  r 


or,.  S  (x. )  = 
rl  cn 


21 


o  r —  r — *  op  p 

r;  (ki,p  ik2>  - pk  j 


kl  k2 


The  Central  i a  \ t  Thoorom  enables  us  to  assume  that 


(63) 


BP(kx  /P^(ic2  )  ^  p2(k1)p2(k2)  +  2  s(k1)p(k2) 


(61*) 


Hence 


Tf, 

2  * — 1 


cn 


2  )_ 


kl  k2 


„  rf 

-2  A-, 


=  2Tfs  p2(k)  +  2 2_ pTk^7Fnr2T 


kl^k2 


rf 

-2  3 


-  2Tfs  P2(k)  ♦  ^(Tf,.  -  k>  P(o)P(k) 


k=l 


or 


cn 


=  2Tf 


,  rr 

-2  , — ,s 


p2(k)  ^  (1  -  ~)  p(o  )p(k  ) 


-2 


k=l 


(65) 


Keeping  in  mind  that  rX/c  is  constrained  to  be  an  integer  we  obtain  that 
.2  0 
"cn  r.j  ( Li-r)‘~ 

2Tf  3  G 


2  -1  I  "|2  rXN  (rt-2r)2 

r  —  I  +  2(1  -  — )  - 


iin  XT  T 

n  N+I 


p  i  T 

1 - sin  rr— r 

n  N+I 
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Tf 

- 

lr=1 


i’f 


“ 

F.-r 

.  -1 

sin 

n 

_ 

Np  (ko)+  Ip. (ko)l 


i  ■  •  '  i '  - '  .-II  , ,  > 

-  -  sin  rr-r  P,  (ko) 


.'I  +  I 


N+I 


,-2r 

•  -1| 
sin  | 

lp^(ko-rX)\  ^ 

|  i'Pn(ko-tX)  +  Ip^(ka-rX)| 

<2n 

x  +  i  rjln 

1  M  +  I  j 

rX  v  . 


(66) 


If  we  now  /:.:.ke  ’.ho  slightly  .rf"ongcr  assumption  that  p(-^-)  is  quite  small, 
the  cross  term  in  the  last  square  can  be  ignored  and  equation  6 L  becomes 
2 


C.n  rj 


ITf 


i  2 

„  2  1 

n  .  -j  I  | 

o  -  sm  cr-y 

2  1  + 1 

1  *  |  (1 

_  rX)  (h“2r) 

T  ’  v  M-  r  ' 

rf 

-.s 


+  2 


!Z  (  1  -  rr 


k=l 


.  -1 
3in 


Npn(kn  )+Ipi(ko)‘| 


M  +  I 


-  sin 


N+I  pi(kc)j) 


i’I 


rX 


■  i  — 2r 
1  vh  -  r  ' 


»  s  o  ,  rX 
2c-'  k+ — 


(1  - 


Tf 


! 


k=l 


.  -1 
sin 

Mpn(ko)+Ici(ko)‘ 

MTI  Pi^k^. 

N  +  I 

23  _L  i  1 

L  J 

\  2 
ih ) 


For  large  decision  times,  this  becomes 

D  2  =  2ff  i^ZLl  (i  _  £  sin-1  — [l  +  \  (1  -  E2l)C-iri£) 
cn  s  Li  n  M  +  I '  [  2  T  i-.  -r’ 


oo 


k-1  5 


where 

1/2  J  -1 
•'cn  (ko)  =<>din 


Mf»n(ko)  +  Ipi(ko  )i 


-sin  ^ 


N  +  I 

Inserting  .uuations  6?  and  31  into  Equation  11  results  in 


f  I  .  |l  /.n  .  -1  I 

!—  Fj_(  ko  /i  f  r  -  -.in 

JJ'  L 


(6o) 


(69) 


_  . 

c :  l  7i  N 


U=F 

,  ,  I /M+I  /  I  \2\  k  N 

h--r-l  +  i  (nTI)  jL  (1-  i— ) 


Pi2U:X) 


Vl-(I/N+I  ) 


k=i'  h-rVl-  ( I/N+I  )^pT^(kX) 


00 


Ui)(l  -  |sin“1nTT¥+?(iTT?)  (1  -  T¥1+C(A  -  tT‘ )acn(ko) 


k=l 


(70) 
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This  last  result  differs  from  th  t  obtained  under  the  assumption  of 
independent  samples  (Equation  jfl  )  ir;  only  two  ways.  First  there  is  the  presence 
of  the  additive  term  in  the  numerator  and  multiplicative  term  in  the  denominator 
as  a  result  of  the  sample  dependency.  The  additional  term  in  the  numerator 
vanishes  for  interference  sources  that  are  far  removed  in  bearing  from  the 
target.  In  fact  it  can  be  shown  by  putting  in  numbers  that  tnis  term  is 
always  negligible,  and  hence  it  will  be  ignored.  The  second  difference  is 
that  there  is  a  tern  missing  in  the  denominator  as  a  consequence  of  the 

ass  .mptioi  that  the  test  statistic  is  normal.  This  term,  which  we  mentioned 

previously,  can  modify  the  SNR  by  as  much  as  1  db  for  I/N  =  10^  depending  on 
whether  or  not  it  gets  multiplied  cy  the  last  term  in  the  denor.unator  of 
Equation  70.  Since  the  cost  of  clipping  for  this  high  I/N  ratio  will  turn  out 

to  be  significantly  larger  than  1  db,  this  possible  error  source  is  not 

particularly  significant  In  fact,  if  we  define  the  signal-to-.i dse  ratio 
foe  uncorrelated  samples  (Equation  38)  as  SHR  ,  then  we  will  taxe  as  a 


possibly  pessimistic  bound 

-i 


SNR  =  SNR 
cn 


cn 


(71) 


oo 

1  *  2Z.(1  -  if  An 

k-1  s 

Correspondingly,  an  examination  of  Equation  21  leads  us  to  conclude 
tha  t  for  t,he  analog  nulling  detector  D  ,  we  can  write 

SNR  =  SNR'  - - - 

an  an 


(72) 


oo 


1  +  2^  (1  -  ^|-)  Pn2(ka) 


k=l 

It  follows  that  th'~  cost  of  clipping  is  equal  to  that  given  in  Table  1  for 
slew  sampling  rates  modified  by  the  improvement  ratio  (I )  which  is 
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I  = 


1  +  2 


I>- 


k=l 


CD 

1  *  2Z(i 

k=l 


fT>  “cn^ 

s 


(73) 


where  .  is  iven  m  .quation  69.  iht  inherent  co^v  of  clipring  is  outained 
by  consid  ring  infinite  ;ar.pling  sates  or 

CD 


o 

I"  = 


I 


(l  -  l  /T  )  p  VodM 


n 


(7U) 


w 

/ 


(1  -  rf/T)  P.cn(H  )dM 


We  can  readily  determine  from  Equation  69  that 


R  (ko  )=rp  (kc  )  - - = - ^ - » - j—  , 

cn  n  n  /,  I  wT  2  .-1  I  , 

U+m)0  -  7sln  n 


(  7?  ) 


n 


f!+I 


1 


i.f  it  is  assumed  that  p.(ka)  =  p  (ko).  The  equality  is  satisfied  for  p  (ko)  =  0 

in  !• 

but  is  very  nearly  satisfied  for  small  values  of  Pn(ko).  The  inequality  is 
a  consequence  of  the  fact  that  hard-limiting  always  costs  something  regardless 
of  the  spectral  shape  of  the  inputs.  However  the  fact  that  both  sides  of 
Equation  7b  are  nearly  equal  for  an  appreciable  range  of  values  of  Pn(ko) 
is  indicative  of  the  fact  that  a  significant  portion  of  the  cost  of  clipping 
for  slow  sampling  rates  is  a  consequence  of  the  sampling  rate  rather  than 
the  hard-limiting. 

k  > 

hssuning  large  decision  times  (1  -  ^-)S1),  and  a  noise  spectra  that  is 

3 

flat  out  to  f  =  wq/2h  CPS  and  zero  elsewhere,  Equation  73  is  plotted  in  Figure 
2  as  a  i unction  of  sampling  rate  with  I/M  as  the  parameter.  The  limiting  case 
(Equation  ~h)  is  apparent  from  this  figure  and  is  tabulated  on  the  following 
page. 
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i/i: 

0 

1 

10 

100 

1000 

I  ( db  ) 

1 

1.5 

2.9 

5.0 

7.U 

in  bio  Li 


h’ote  also  the  "improvement"  of  the  PGA  relative  to  the  analog  devices  corresponds 

to  the  1  db  given  in  Table  L  for  I/!’  =  0.  Thus  tnis  table  also  represents  the 

imrrover.ent  of  D  relative  to  the  PGA  detector  if  all  the  values  are  reduced 
cn 

by  1  db.  ,-.s  an  example  let  us  assume  an  array  of  UO  hydrophones.  '’aole  5 

gives  the  net  irr.orovement  of  IJ  relative  to  the  P  :«  as  well  as  the  net  cost 

cn 

of  clipping  r*  lative  to  0_  . 


r/r: 

0 

1 

10 

— 

100 

— 

1000 

J  (  db ) 

-1 

it.  2 

12.5 

20.6 

28.8 

C  ( db  ) 

1 

1.7 

3.7 

6.2 

8.2 

Table  5  (M  =  UO) 

Equation  73  is  also  plotted  in  Figure  3  for  a  noise  spectrum  that  is  not 
pre-whitened  but  falls  off  at  a  rate  of  6  db/octave  beyond  f^/Li  and  is  cut-off 
sharply  but  not  ineally  at  f^cps .  While  the  inherent  cost  of  clipping  is  most 
aonrooriately  determined  for  pre-whitened  spectra,  Figur’d  3  indicates  that  the 
results  of  'able  d  are  relatively  insensitive  to  the  spectral  shape. 

Finally  we  are  concerned  with  the  cost  of  clipping  when  sampling  at  the 
rather  arbitrarily  determined  "minimum  reasonable"  sampling  rate  f  .  For 
I/H  -SF  1,  i  yquist  rate  sampling  is  adequate  and  there  is  no  improvement  for 
"rapid  sampling".  for  lar^e  values  of  I/N,  the  improvement  is  approximately 
1  db  less  than  that  given  in  "able  U.  These  results,  obtained  from  Figure  2, 
are  oased  on  the  assumption  that  the  array  can  be  aimed  directly  at  the  inter¬ 
ference.  When  thr  interference  is  midway  between  permi ssable  bearing  angles, 
the  sampling  cost  ie  of  the  order  3  db  for  D  and  li  db  for  D  .  Hence  for 
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this  \itor  case,  the  total  cost  of  sampling  and  clipping  is,  if  anything,  less 
than  given  in  Table  «lso  for  otner  ti.an  ideally  flattened  spectral  shapes, 
the  improvement  in  the  cost  of  clipping  lor  sampling  at  'he  f  rate  is  greater 
ti.an  that  obtained  from  figure  2.  For  example  the  spectral  shape  that  determines 
figure  3  gives  an  improvement  factor  at  the  f  sampling  rate  that  is  almost 
identical  to  the  results  of  Table  U,  even  when  thi  array  is  aimed  directly 
at  the  int;  rference.  A  reasonable  rule  of  thumb  for  the  cost  of  clipping 
for  moderately  large  arrays  and  for  sampling  rates  f  >  f  flight  be 
T<  e  cost  of  clipping  is  a;- proximatcly  2  db  for  an 
interference  to  ambient  noise  power  ratio  I/N  »f 
unity  and  increases  almost  linearly  with  I/N  at 
a  rate  of  2  db  per  decade. 

VI T .  On  the  Foint  .Source  Assumption  and  Hydrophone  Tolerances 

The  sampling  rate  discussion  of  Section  V.  gives  us  a  convenient  frame 
of  reference  with  which  to  decide  whether  or  not  a  source  can  ue  regarded  as 
a  point  source.  Let  us  assume  that  the  detector  is  sampled  at  the  rates  given 
in  Table  3  which  are  just  adequate  to  null  out  a  point  source  with  a  given 
I/N  ratio.  The  resulting  bearing  quantization  determines  that  maximum  length 
that  a  target,  at  a  particular  range  and  with  a  given  I/N  ratio,  can  have  and 
still  ie  regarded  as  a  point  source.  Thus,  if  the  length  of  the  target  is  no 
greater  than  the  product  of  the  range  with  the  minimum  change  in  bearing 
angle,  then  there  is  no  ambiguity  as  to  its  direction  and  it  is  for  all 
practical  purposes  a  point  source.  Another  way  of  looking  at  this  is  to  realize 
that  such  a  target  can  in  principle  bu  nulled  out  at  a  cost  that  is  comparable 
to  the  avv.ra,.e  cost  of  sampling  at  the  f  rate  when  nulling  out  a  perfect 
point  source. 
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If  the  interference  is  broadside,  then  o  =  1/f  =  d/c  sinQ.  Hence 

s 


AO  =  c/f  d,  and  therefore 
'  sn  ’ 


f 

^  c  r  o 

q*  f  d“ 
o  sm 


(76) 


where  q  is  the  effective  length  of  the  target  and  f  /f  is  given  in  (’able  3 

s  o 

as  a  function  of  I/K.  For  endfire  interference  Equation  76  becomes 


,  2  f 

1  c  r  o 

2  f  d  T 


1 

2 


(77) 


|_  o  sm  J 

Let  us  assume  that  the  inputs  are  pre-whitened  out  to  5  kc  (  j  =  1  )  and 

o 

a  hydropnone  spacing  of  2  feet.  The  maximum  effective  target  length  is  given 
in  Table  6  as  a 


Broadside  < 


Endfire 


Thus  if  the  target  dimensions  are  the  effective  propellor  dimensions,  it  can 
come  as  close  as  1000  and  supply  an  I/N  ratio  of  30  db  and  still  be  re;  arded 
as  a  point  source.  In  fact  for  a  20  db  I/N  ratio  an  entire  submarine  at 
1000  yards  can  be  regarded  as  a  point  source.  We  conclude  that  the  point 
source  assumption  for  interfering  shipping  is  not  unreasonable. 

Finally  let  us  consider  the  allowable  tolerance  on  the  hydrophone  spacings. 
in  order  to  null  out  an  endfire  interference  the  maximum  delay  per  hydrophone 
of  d/c  seconds  is  required.  If  the  actual  hydrophone  spacinge  is  Ad  feet 
different  from  the  assumed  spacings,  then  the  required  delay  would  be  in  error 
by  Ad/c  seconds.  If  we  now  require  that  this  error  be  no  greater  than  l/2f 

sm 


function  of  range  and  I/N  ratio. 


"'vMi 

r 

1 

10 

100 

1000 

2000 ' 

555.' 

175.5' 

55.5' 

17. 6* 

1000 ' 

277. 5* 

87.8' 

27.8' 

OO 

• 

CO 

2000 ' 

7U5. ' 

109.  ' 

235.5' 

132.5' 

1000 ' 

372.$' 

209.5' 

118.  ' 

66.  ' 

Table  6 
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seconds,  the  nulling  will  still  e\>  effective  even  for  endfire  interference. 
Thus  a  reasonably  tight  bound  on  the  hydroonono  spacings  might  be 


f 

lids  6  -p  inches.  (?8) 

o  Ism 

Hence  if  we  assume  that  the  inputs  are  processed  out  to  5  kc,  the  tolerances 
in  the  hydrophone  spacings  needed  to  null  out  an  endfire  interference  whose 
power  is  I/N  that  of  the  ambient  noise  is  given  by 


I/to 

r  - 

1 

10 

100 

1000 

Ad 

3" 

1'' 

.3" 

.1" 

Table  7 
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Surnrna  ry 

The  likelihood  ratio  test  for  detecting  a  weak  signal  in  the  presence 
of  a  strong  point  source  interference  has  been  evaluated,  given  that  the 
inputs  have  been  reduced  by  hard-limit ing .  An  implementation  which  very 
closely  approximates  this  test  and  is  in  fact  asymptotically  optimum  as  the 
interference  to  ambient  noise  power  ratio  increases  is  presented.  This 
procedure  aims  the  array  in  the  direction  of  the.  interference  and  subtracts 
the  average  value  of  the  inputs  from  each  of  the  inputs  before  proceeding  with 
the  standard  DIMUS  processing.  As  this  procedure  is  asy:nptoticnlly  optimum 
for  analog  inputs  as  well,  it  is  conjectured  to  be  optimum  for  arbitrary 
quantization.  It  is  shown  that  the  degradation  due  to  hard-limiting  indicated 
in  report  26  is  indeed  real  and  that  a  relatively  large  number  of  quantization 
levels  are  needed  to  significantly  reduce  this  cost  for  large  interf erence- 
to-noise-racios . 
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I  ■ _ Introduction 

This  report  deals  with  the  detection  of  a  sonar  target  in  the  presence 
of  a  strong  point  source  interference  with  the  constraint  that  the  inputs 
are  reduced  by  arbitrary  quantization.  It  will  be  assumed  throughout  that 
the  hydrophone  inputs  in  the  absence  of  the  signal  and  interference  targets 
are  statistically  independent.  This  problem  was  considered  without  the  data 
reduction  in  reports  17  and  21.  The  likelihood  ratio  test  was  evaluated  and 
it  was  shown  that  the  interference  was  effectively  removed.  Furthermore, 
a  suboptimum  implementation  called  the  nulling  detector  was  proposed  which 
closely  approximated  the  likelihood  ratio  test.  The  nulling  implementation 
was  further  analyzed  in  report  26  under  the  assumption  chat  the  inputs  were 
hard-limited.  While  this  ins trument a 1 1  op.  represents  a  vast  improvement  over 
the  standard  polarity  coincidence  array,  the  results  indicated  that  the 
degradation  in  detectability  due  to  hard-limiting  is  not  negligible  for  strong 
Interfering,  targets.  r.i  order  to  determine  whether  or  not  this  degradation 
is  inherent  in  the  problem,  one  must  determine  the  performance  of  the  like¬ 
lihood  ratio  test  for  the  clipped  data.  If  the  degradation  is  indeed 
significant,  it  will  be  useful  to  estimate  the  improvement  in  detectability 
for  multilevel  quantization.  This  result  should,  In  addition  to  indicating 
the  advisability  of  going  to  multilevel  quantization,  indicate  the  accuracy 
required  for  the  detector  employing  analog  data. 

There  is  no  known  approach  for  the  general  solution  of  the  likelihood 
ratio  for  quantized  inputs.  The  only  known  procedure  is  to  assume  that  the 
inputs  are  sampled  slowly  enough  so  that  the  samples  are  independent  (Nyquist- 
rate  sampling  for  gaussian  inputs)  .  Faster  sampling  rates  are  required 
however  to  achieve  a  low  cost  of  clipping  and  to  permit  reasonable  accuracy 
in  bearing  measurements.  In  fact,  it  was  indicated  in  report  26  that  unless 
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cne  sampled  rapidly  enough  Che  major  source  of  degradation  in  the  nulling 
detector  is  due  to  the  sampling  rate  rather  than  the  quantization.  Once  the 
optimum  detector  for  independent  samples  has  been  found,  however,  it  can  be 
analyzed  for  rapid  sampling  rates.  There  are  a  number  of  reasons  whv  this 
procedure  is  not  thought  to  be  restrictive.  In  the  first  place  there  is  no 
a  priori  reason  to  assume  that  the  form  of  the  optimum  detector  would  be  a 
function  of  the  sampling  rate.  For  example  in  the  case  of  analog  data,  the 
only  difference  is  the  presence  of  pre-whitening  filters  at  the  inputs. 
Secondly  even  if  the  optimum  deLectoi  did  depend  on  the  sampling  rate,  one 
would  need  precise  knowledge  of  the  ambient  noise  power  spectrum  in  order  to 
implement  it.  The  complexity  of  such,  a  system  would  rule  it  out  even  if  such 
knowledge  were  available  and  it  generally  is  not.  V.'e  will  therefore  find 
the  likelihood  ratio  for  independent,  hard-limited  samples  of  the  hydrophone 
inputs.  It  will  be  assumed  that  the  array  can  be  aimed  either  directly  at 
the  interference  or  the  target  even  though  the  sampling  rate  required  by  this 
assumption,  as  discussed  in  report  26,  is  inconsistent  with  the  Independent 
sample  assumption.  Thus  the  detector  wc  Te  seeking  is  to  be  optimum  for 
independent  samples,  but  it  will  in  fact  be  operated  with  dependent  samples. 
II .  _  T h e  L i k e llhcod  Ra_ti o_ f  or_  We ak_  S ig nals 

A  graphical  representation  of  all  of  the  available  data  is  shown  in  Fig. 
1  where  it  is  assumed  that  the  array  is  aimed  directly  at  the  target.  The 
assumption  that  the  array* could  be  aimed  at  the  interference  is  indie. ted 
by  the  fact  thaf^  the  interference  wave  front  intersects  M  data  points 
simultaneously.  It  follows  from  the  assumptions  that  the  noise  inputs  are 
all  statistically  independent  and  the  samples  in  any  channel  are  independent 
that  the  correlation  between  samples  along  any  vertical  lipe  (Yixed  i,nstant, 
of  time^  is  due  solely  to  the  common  additive  signal.  Given  that  the  data  is 
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ha  rd- limited ,  the  information  available  at:  to  the  presence  of  the  signal 
is  contained  in  the  number  of  samples  along  a  vertical  axis  that  have  the 
same  sign  where  all  of  the  samples  are  suitably  weighted"  by  the  number  of 
samples  that  have  the  same  sign  along  the  diagonal  axis  which  contains  the 
common  additive  interference.  l.'e  will  assume  that  the  number  of  samples  N 
is  sufficiently  large  Jn  .->  M(M-1)1  so  that  the  "corners1  of  Fig.  1  where 
the  diagonal  lines  have  fewer  than  M  elements  can  be  ignored.  The  optimum 
implementation  will  involve  aiming  first  at  the  interference  and  then  after 
some  processing,,  aiming  at  the  target.  Hence,  the  data  array  will  really 
be  in  the  form  of  a  rhombus  with  the  corners  missing. 

Let  us  denote  the  probability  that  a  sample,  viith  an  additive  sipnal 
of  value  £,  is  positive  and  x (  of  the  M  samples  along  the  common  inter¬ 
ference  diagonal  are  positive  by  P  +  (x  ,£;).  Hence 


oo 


P+(x1,C) 


f 1 

(m-i  I ' 

J 1 

CD 

,xi-1J[ 

x , - 1 ! 


Fn+s( 


-u> 


M-x , 


1-F  (-U-OI  Mu) 
n  i 


du 


(!' 


where  F(  )  is  the  gauasian  cumulative  distribution  function  and  f^(  )  is 
the  density  function  of  the  interference.  Similarly  P  (x^f)  is  the 
probability  that  the  sample  is  negative  and  x^  of  tbo  sa-ples  along  the 


interference  diagonal  are  positive- 

,  f 

(-u)l 


P  (*l.O  = 


oo 

f  I M-l' 


J  lXi 

-oo 


i-F 


n+S 


F  ,  (-u) 
n+s 


1  M-x^-i 


F  (-u-r)  fT(u)  du 
n  v  t 


(2) 


All  that  matters  is  the  number  of  samples  along  a  vertical  axis  that  are 
positive,and  the  ordering  of  them  is  irrelevant.  The  probability  that  r 
samples  are  positive  with  an  arbitrary  distribution  along  the  interference 
diagonals  is  given  by 
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oo 


-oo 


(3) 


where  the  subscript  K  refers  to  the  alternative  of  signal  present.  Under  the 
hypothesis  of  no  signal 


M-r 


H(r-x)  =  r  TT  P  <XH  .°)  IT  P  °) 

i=l  j  =  l  ] 


The  likelihood  ratio  L(r,x)  is  given  by 


Pk  (r  >2i) 

L(r,x)  - 


Pj.(r  ,x) 


(A) 


(5) 


and  depends  on  the  signal-co-noise  ratio  S/N  among  other  things. 

The  standard  procedure  is  to  consider  a  locally  optimum  detector  that 
is  optimum  for  vanishingly  small  signa.l-to-noise  ratios.  Hence  we  in  fact 
implement  the  fi.:st  two  terms  only  of  the  expansion 

T  •  ( r  ,x  ,S/N)  =  L(r  ,x  ,0)  +  I  ’  S/N  + 

1  ;  S/h'-O 

We  will  denote  these  first  two  terms  by  TT(r,x).  This  likelihood  ratio  will 

be  evaluated  by  expanding  the  probabilities  P+(x^,£)  and  P  (x^.^i  and  retaining 

2 

only  those  terms  of  orkr  (S/N)  or  £  and  then  repeating  this  procedure  for 
P  (r  x) .  It  may  be  that  when  we  are  done,  some  of  the  terms  of  L(r,x) 

In.  ' - 

will  not  represent  cross-correlation  information  and  hence  are  useful  only 
when  the  noise  power  is  a  priori  knowledge.  Ignoring  such  terms  would  in 
ft  ct  represent  no  real  loss  of  information. 

Expanding  and  then  collecting  all  of  the  terms  in  Eq.  (1)  is  straight¬ 
forward.  Retaining  only  the  pertinent  terms,  one  obtains: 

P+(x±,0  ^>P(xi,0)  +  P*^)  S/N  +  PB+(x.,)E  +  Pc+(xi)E2  ,  (7) 


where 
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p<v0>  ■  ft!" 


x,,  ,  fl-x 


l-Fn(-u)]  1[pn(-u)]  1  fj(u)  du 


(8) 


PA+0<1)  -Bi(“JM2-i  /  [l-p„(-“>]  lfn(-U>  "  fI<U)  dU 


-oo 

oo 


'  FT  IxJ  ~2~'  J  [1-Fn(-u)]  [Fn(-u)J  1  fK(_u)  u  fI(u)  du  ’ 

-oo 


(9) 


and 


oo 

+  Xi  'M  1  Ti  1  Xi-1r  -|  M-x 

PB  <*i>  “  FT  Ixl  J  [1~Fn(_u)J  [Fn(-U>J  fn(-u)  fI(u)  du  ■  (10) 

“OO 


Pc+<1‘1>  -r  IxJ  /i1-F„<-“>]Xl  1[Pn<-'‘>]H,‘1lN  V'“>  •  !ll) 


Similarly ■ 


M“X 

P“(x1,C)-»  [p(xt,0)  +  pa“(x,)s/n  -  PB”(x±)c  -  Pc~(x.)(;2  |  , 


(12) 


where 


x, ...  M-x  -1 
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OO 
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M-x . -1 


-  sr  !”J  j  ‘[V-’l  1  sV->  'i«  ■  <15> 


If  one  now  substitutes  Eq  .  (7)  into  the  first  product  term  of  Eq .  (3)  and 
then  retains  nn,y  che  pertinent  terms,  one  obtains 

r  +  r  -T  r 

FT  ?  <x  .  0  — >TT  P(x.,0)  +  )  P  (x  )  TT  P(x,,0)  S/N 

i=l  1=1  1  ,  ,  id-1  3 


i  =  l 


r  r 

—  -*  r  _ ,  r 

'  pnTx  )  Ti  P(x  ,0X  +  )  P  +(x  )  FT  P(x  0)  C 

“  J/l  J  ^  3*1  J 


i=l 


+  2  PB+<*J>  Tf  P<«k-« 


or 


,  r  •>  ,  ,  r_,+  . 

r  r  r-  p.  \x  )  (x  ) 

FJ^P  (xi/,)^TT  ?  (X1 1°)<  1  +  S/N  ^  pTx~0)  'K2^P(x1,0) 

'  i=l  1=1 


+  f  2  v  c  x*  +  l i  »  b  (xi} 

■  L  p(x,,0)  +  2  2_-/_.P(x  ,0)  P  (x .  ,0) 

1=1  1  idj  J 


Similarly 

M-r 


M-r  M-x 


M  r  pT(xJ  M_r 


(16) 


(17) 


pn  (x±) 

of 


TT  My)  -  TT  ^  '(V»>  1  *  L  ^  re- 

J  J  3  {  j=l  3  j=l  1 


,-r  yy>  f.2vy-y<*t>  p,'(»,) 

£—>  P(x  ,0)  2  Z_  Z_P(x. ,0)  P(x  ,0) 

j=l  T  if}  3 

It  follows  from  Eq .  (A)  as  well  as  Eqs.  (17)  and  (18)  that 

r  M-r  M-x 

pH(r,x)  -TT  TT  P(^.0)  TT  — ^  P(x  0) 

H  “  i=l  j=l  xj  T 


(18) 


(19) 
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Combining  Eqs.  (3,5,17,18,  and  19)  one  obtains  for  the  likelihood  ratio 


L(r,x)  =  1  +  S/N 


r  P.V) 


M-r 


V  rA  +  r  pa"(xi)  V  Isjv 

L  P(X  ,0)  L  P(x.,0)  "  NZ-  L  P(x,  ,0)  P(x.  .0) 

i'l  i=l  1  i=l  1=1  1  3 


r  M'r  ?_  (x.)  P  "(x,) 
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r  r  P„+(xJ  Pn+(xJ 


+  Nr  c  ^±'  _  r1  c  N^rrB  vAiy  b 
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i=l  1=1  1  i*j  1  3 


M-r  M-r 

fl  I 

if) 


PB  (xi3  PB  ^X13 
P(x1,0)  P(xj(0) 


(20) 


It  is  now  conjectured  that  the  terms  involving  P  (x  ) ,  P  (x  ) , 

Pc+(xi) ,  and  P^,  (x^)  either  vanish  or  cancel  themselves  out  depending  on 
the  value  of  the  interf erence-to-noise  ratio  (I/N) .  It  Is  not,  however, 
necessary  to  prove  this  conjecture,  since  we  can  argue  that  these  terms 
should  be  ejected  in  the  event  that  they  do  not  cancel  themselves  out. 

One  notes,  by  tracing  their  origins,  that  these  relatively  few  terms  (for 
large  arrays)  do  not  represent  cross-correlation  information  and  hence  can 
only  be  useful  when  the  noise  power  is  a  priori  knowledge.  They  are  completely 

of  a  total  of  M'J  of  the  standard  array  detector 
that  measure  the  power  of  the  hydrophone  inputs  rather  than  the  correlation 
between  them.  When  the  noise  power  is  not  known,  these  terms  can  onl,  ad 
to  the  variance  of  the  test  statistic  and  hence  they  should  be  omitted 
(see  reports  18,22,  and  23).  Since  one  indeed  never  has  a  priori  knowledge 
of  the  noise  power,  the  useful  likelihcod  ratio  is  given  by: 

{r  r  M-r  M-r  r  M-r 

Y_  NV(x1)P+(xj)  +  Y_  Z  P‘(xi)P"(xj)-2^  ^  P+(x1)P‘(xJ| 
itH  i?H  1  J  (2D  I 


where 
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+  a  pn+(x 
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)  (22) 


P'(Xi>  "  ?fx. ,0)  ■  “oc 
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fj(u)  du 


(23) 


1’  rr  X  ..M-X 
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Let  us  now  evaluate  the  likelihood  ratio  for  the  case  of  no  Interference 
'I-^?0).  Substituting  for  f^(u)  in  Eqs.  (22  and  23)  the  delta  function 
6 (u) ,  one  obtains  by  inspection  that 

P+(Xl)  -  P'(Xi)  =  2  fn(-0)  ^  .  (24) 

It  follows  from  Eq.  (21)  that 

L(r,x)  1  — >  (l  +  ~  |-  [(r2-r)  +  (M-r)2  -  (M-r)  -  2r(M-r)j| 

-  j* +  i |[<2r"a>2  -  M]j  •  <“> 


This  is  exactly  the  result  obtained  in  report  6  for  the  standard  detection 
problem,  and  it  therefore  acts  as  a  partial  check  on  Eq.  (21).  It  is  known 
that  the  DIMUS  array  implements  the  likelihood  ratio  of  Eq.  (25)  See  report 
22  for  example  . 

III.  An Implementation  of  the  Likelihood  Ratio 

It  is  informative  to  interpret  the  meaning  cf  the  weights"  given  by 
Eqs.  (22  and  23).  In  the  case  of  no  interference,  we  note  that 

P+  ■  Prob  /  n  tS  0/n  >  0 
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and 


P  (x^)  “  Prob  {  n  c  0/n  <  0  (26) 

These  weights  are  obviously  pertinent  for  the  weak  signal  case,  since  they 
are  the  probability  that  a  vanishingly  small  signal  determines  the  sign  of 
the  samples,  ifhen  interference  is  present,  we  recognize  that 

[l  -  F  (-u)j  1[f  (-u)l  1  fT(u) 

-  =  Prob  i  i(t)  “  u/x  ^  .  (2/) 


M-x , 


[l  -  Fn(_u)]  i[^Fn(--u)|  1  fj(u)du 


It  follows  that 


P  (x^)  **  Prob 


P  (x^)  =  Prob 


cf-  -i (t) /n  >  -i(t) ,x. 


;  -i (t) /n  <  -i ( t)  ,x 


(26) 


Once  again  we  have  the  probability  that  a  vanishingly  small  signal  determines 
the  sign  of  the  samples.  Note  that  if  we  used  a  standard  DIMUS  system  that 
did  not  use  the  information  contained  in  the  vector  x..  the  weights  would  be 


^2/^n+i)'  which  are  not  the  optimum  weights  given  by  Eqs.  (22  and  23). 

It  is  intuitively  clear  from  the  form  of  the  likelihood  ratio  ^ Eq .  ( 2 1 ) j  , 
however,  that  the  optimum  implementation  is  a  DIMUS  system  that  Is  preceeded 
by  some  type  of  pre-processing  which  is  based  on  the  number  of  positive  samples 
(x^)  along  the  common  interference  fronts.  Insight  is  gained  into  the  type 
of  pre-processing  by  considering  some  special  cases. 

It  is  eauily  shown  from  Eqs.  (22  and  23)  that  when  x^  =  M,  P+(x^)  — — y 

0  as  I/N — ;>oo  and  when  x^  *»  0,  P  (x^) - * — J>0  as  I/N — >oo  .  It  follows  that 

when  all  of  the  samples  along  an  interference  fr^rt  have  the  same  sign,  they 
contribute  very  little  information  about  the  presence  'of  the  signal  and  for 
large  Interference  sources,  they  should  be  ejected.  It  can  be  numerically 

demonstrated  that  the  optimum  '‘weights"  become  negligibly  small  when  the 
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interference- Co-noise  ratio  I/N  is  large  enough  so  that  the  probability  of 
the  events  x^  =  M  and  x^  =  0  become  significant.  Thus,  it  costs  very  little 
to  eject  such  samples  regardless  of  the  interference-to-noise  ratio.  Recall 
that  in  the  analog  case  (Report  17),  the  asymptotically  optimum  (I/N  -yoo) 
implementation  very  nearly  approximates  the  optimum  for  large  arrays  ^intro¬ 
duces  a  cost  of  one  .'.ydrophone  for  I/N  =  0  j  and  does  not  require  a  priori 
or  estimated  knowledge  of  the  interference-to-noise  ratio.  It  will  be  oeen 
that  this  is  also  the  case  when  the  inputs  are  hard-limited.  Ic  can  also 
be  numerically  demonstrated  that  for  large  values  of  I/N,  the  probability 
of  the  event  x  =  M  or  0  is  i-t,  where  t  is  of  the  order  /n/I !  The  removal 
of  such  a  .large  percentage  of  the  samples  constitutes  the  significant 
degradation  due  to  hard-limiting  that  was  observed  in  Report  26. 

It  is  also  informative  to  consider  the  average  weight  along  a  common 
^nteference  front.  This  average  weight  w (x^)  is  calculated  from  the  formula 

W(x1>  -  ^  [x^*^)  ~  (M  -  *±)P  (x^J  .  (29) 


It  turns  out  that  W(x^)  — c?0  as  I/N— >>oofor  all  values  of  x^.  This  is  difficult 
to  show  rigorously,  but  it  can  easily  be  shown  to  be  approximately  true  by 
considering  a  suitable  approximation  to  P+(x^)  and  P  (x^)  .  This  result  has 


already  be».;n  demonstrated  for  x^  -  M  or  0.,  For  other  values  of  x^,  it  can 
be  shown  that  in  the  limit  as  I/N — goo,  the  probability  of  the  event  fi(t)  ° 


u/x^l  as  given  by  Eq .  (27) 


is  an  even  function  about  its  maximum  which 


occurs  for  that  value  of  u  such  that 


(30) 


f  (-u)  f  n  (~u) 

The  functions  r— — ; — and  rr—, - r-  are  very  nearly  linear  functions  and 

[1-7n("u)]  Fn("u) 

hence  P+(xi)  and  P~(xi)  can  be  jpproximated  by  these  functions  evalusted  at 
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u  ■*  u  (see  Eqs.  22  and  23).  Thus 
max  M 


^  2^  eXpf  2  [£  (  M~) j 


’(xi}  ~:[i  Z(M^x^)exp ;  \  5"  '(m 


For  these  approximate  weights,. Eq.  (29)  Is  identically  zero.  The  exact  weight, 
differ  in  the  same  direction  so  that  U(x^)  is  still  identically  zero.  It 


is  interesting  to  note  that  u  is  the  maximum  likelihood  estimator  of  the 

max 


interference  and  also  that  P  (x^)  ar.d  P  vx^)  for  ■=  M/2  ere  what  they  should 
be  for  I/N  =  0. 

The  previous  results  suggests  an  implementation  suggested  by  Anderson 
for  the  case  of  analog  inputs  J  which  consists  of  steering  the  array  .at  the 
interference  and  then  subtracting  the  average  value  of  the  samples  from  each 
of  the  samples.  This  implementation  would  weight  the  samples  in  such  a  way 
that  the  average  weight  would  be  indeit  ically  zero  and  in  the  event  x^  -  M 
or  0  the  weight  would  be  zero.  If  x,  samples  are  positive,  the  average  valve 


would  be  1/M|x^-(M-x4)  =•  2x^/M  -1.  It  follows  that  the  modified  sarnies 

would  have  the  values  2(l-x^/M)  for  positive  samples  and  -2w  /M  for  negative 
samples.  This  procedure  therefore  introduces  multiplicative  weights 
2(l-xi/M)  and  2x  /M.  Regardless  of  the  value  of  I/N.  the  expected  value  of 
x^  is  equal  to  M/2  (provided  only  that  t  =  0)  and  hence  the  expected  value  of 
the  weights  is  equal  to  unity.  This  modified  DIMUS  system  is  analyzed  in 
Appendix  A  where  if.  is  shown  that  for  weak  input  signals,  the  output  signal- 
to-noise  ratio  becomes 


2,  _S 
ti|  2tt  N+ 


_  ( 1-1 /M)_ 

t  r.  2  V  -l  "i  i.r,~  rlr 

L1  rr  sin  N4I.V1  ’  3M 
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Note  that  for  large  arrays,  the  cost  of  this  implementation  is  quite  small 
(about  4/3  hydrophone)  for  vanishingly  small  interf erence-to-noisc  ratios 
even  though  this  implementation  is  conjectured  to  be  optimum  only  asymptotically 
as  7. /N  — >oo  .  For  large  values  of  I/N  we  can  replace  sin  ^(I/N+I)  by 


u/ 2-  (/2N/N+I  (see  report  2u)  .  Hence  Eq .  (32)  become'' 


SNR 


I/N  oo 


rw  1 


Hi/ 


Jl-l/H) 


(n+i)  n  x.' 

1/ 1  3M  2 

V  3M 


(33) 


A  comparison  of  this  result  with  report  2b  shows  that  the  nulling  implementation 
approximates  this  asymptotically  optimum  implementation  in  the  identical 
fashion  as  in  the  analog  case  fan  additional  cost  of  0.88  db  ;  .  As  a  result 

L  J 

there  is  no  point  in  analyzing  this  detector  for  rapid  sampling  rates  since 
the  analysis  for  the  nulling  detector  carried  out  in  report  26  should  suffice. 

Let  us  now  consider  the  discrepancy  between  the  weights  of  the  proposed 
scheme  and  the  optimum  weights  of  Eqs .  (22  and  23).  Let  us  replace  the 
weights  of  the  proposed  scheme  ^2(l-x^/M)  and  2(x^/M)  j  by  the  weights 
Kf(x^)2(l-x^/M)  and  Kf (x^) 2(l-x^/M) .  These  modified  weights  would  also 
satisfy  the  special  cases  considered.  The  constant  K  would  of  course  not 
effect  the  performance  of  the  detector.  It  is  clear  from  Eq .  (31)  that  these 
modified  weights  would  approximate  the  optimum  weights  of  Eqs.  (22  and  23) 
provided 

/  l  irl  *1  12 


f(x±)  =  — - - -  exp  , 


- 1  [r1  (i/*1 


(34) 


This  function  is  however,  very  nearly  equal  to  unity  for  all  values  of 
x^/M  and  is  plotted  in  Fig.  2.  In  fact  the  difference  between  f(x^)  and  unity 
is  a  result  only  of  the  fact  the  Eq .  (31)  gives  approximate  weights  rather 
that  the  precise  weiphts.  The  pioposed  weights  are  proportional  to  the  optimum 
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weights  of  Ego;  (22  and  23) .  It  follows  that  the  modified  DIMUS  system 
represents  an  asymptotically  optimum  implementation  (I/N — ^oo  )  of  the 
likelihood  ratio. 

IV .  On  the  Improvement  in  Detect abili ty _ f o r  Mult i 1 eve  1  Quan t i z atlon 

The  asymptotically  optimum  implementation  is  the  same  for  analog  data 
£ Anderson's  proposal j  as  well  as  for  hard-limited  data.  It  is  therefore 
reasonable  to  assume  that  this  implementation  is  optimum  for  an  arbitrary 
degree  of  quantization.  ITiile  an  analysis  of  the  general  case  is  quite 
difficult,  a  figure  of  merit  which  gives  an  estimate  cf  the  improvement 
in  detectability  for  multilevel  quantization  is  rather  easily  obtained. 

This  figure  of  merit  seems  to  he  accurate  enough  to  give  reasonable  quanti¬ 
tative  results. 

For  binary  quantization,  the  samplps  fall  into  two  classes  those  which 
have  the  same  sign  as  the  most  common  sign  along  the  interference  front  and 
those  which  are  switches".  The  samples  are  then  weighted  according  to  which 
class  they  belong  to  and  the  best  estimate  of  the  interference.  Assuming 
that  the  interference  signal  has  a  zero  mean,  the  probability  of  a  "switch1' 
is  identically  equal  to  the  probability  that  the  sample  has  the  opposite  sign 
of  the  interference.  Hence  the  probability  of  a  switch  (P  )  is  given  by 

co 

P0  °  Prob  j^n+s+i>0/i<o]  =>  2P  j^n+s>-i ,  i<C)j  =  2  I  ^l-Fn+s(u)j  fj(u)  du  .  (35) 


o 

Assuming  all  the  processes  are  gaussian,  this  probability  becomes 
for  large  interference  to  ambient  noise  power  ratios; 

oo 


P 

s 


I/N- 


i  co 

-i 


if  iM 


|n+s 


du 


I  ,/K+s'  S/N— >  0,  I  IN 
T!  Y  I  *  ry  I 


asymptotically 


(36) 


The  optimum  pre-processing  effectively  reduces  the  number  of  samples  along 
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any  interference  front  to  twice  the  number  of  switches.  That  is  the  weights 
of  the  switches  vary  from  1  to  2  as  the  weights  of  the  non-switches  vary 
from  1  to  0.  One  can  postulate  a  "cost'”  in  terms  of  the  average  "effective 
reduction  in  sample  size".  This  "cost"  becomes  asymptotically 


,T  •  (37) 

n  V  21 

This  cost  is  essentially  equivalent  (to  within  U  db)  of  the  actual  asymptotic 

cost  in  detectability  obtained  from  Eq .  (33)  for  large  arrays. 

This  figure  of  merit  can  be  generalized  for  an  arbitrary  degree  of 

quantization.  The  asymptotic  "cost"  converges  nicely  to  unity  as  the  number 

of  quantization  levels  increase  and  it  is  conjectured  that  the  results 

can  be  used  quantitatively  as  well  as  qualitatively.  Let  us  assume  that 

the  inputs  are  now  quantized  into  2  q  equal  intervale  as  in  the  figure 

below.  Further  assume  that  the  interference  lies  in  one  of  the  closed  intervals 

I  where  (k-1) A< I,,<k&  and  that  the  samples  are  assigned  values  corresponding 
k  K 

to  the  center  of  che  intervals. 


The  samples  fall  into  three  categories-  those  which  "switch"  to  the  left, 
the  switches  to  the  right,  anu  the  non-switches.  The  nulling  procedure  will 
modify  the  weights  of  the  switches  to  -1  and  1  and  remove  the  non-switches. 
Hence  for  such  interference  fronts,  the  optimum  processing  "effectively' 
reduces  the  number  of  samples  to  the  switches.  However,  for  those  inter¬ 
ference  fronts  which  correspond  to  the  two  open  intervals,  the  pre-processing 
effectively  reduces  the  number  of  samples  to  twice  the  switches. 
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sampling  rates.  Reasonable  sampling  rates  reduce  the  cost»  in  db> very  nearly  in 
half .  Hence  doubling  the  no.  of  quantization  levels  reduces  the  cost  from 
approximately  6  db  to  approximately  db.  More  than  eight  3  ;vels  are 
needed  to  reduce  the  cost  to  3  db ,  Uhile  these  results  are  only  rough 
estimates  of  the  actual  cost,  it  is  clear  that  a  relatively  large  no.  of 
quantization  levels  are  needed  to  reduce  the  inherent  cost  of  hard-limiting 
to  some  minimal  value  for  large  interf erence-to-noise  ratios  (  20  db) . 

However,  the  accuracy  required  for  analog  processing  is  not  very  restrictive. 

It  is  only  necessary  to  knov;  the  value  of  the  samples  to  within  one  standard 
deviation  of  the  ambient  noise.  The  only  equipment  problem  is  the  rather 
stringent  linearity  requirement. 
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Appendix  A 

The  proposed  scheme  calculates  the  following  test  statistic* 


Tf  s 


=  V  \ 


M  M 

V| 

U. 

1  i*J 


x1(k) 


-  y, 


vk> 


where 

x1(k1  «=  sgn^u1(t+kT)j 

and 


(A-l) 


M 

y±  -  k  Y  xr(t  +[  k‘»n!(r-i>]T)  > 

r=l 

where  u^(t)  is  the  input  to  the  iCa  hydrophone  at  Lime  t,  f  is  the  sampling 
rate,  and  T  is  the  decision  time.  The  quantity  is  the  average  of  the 
hard-limited  samples  along  the  interference  front  associated  with  u^Ct+kr  ). 
It  has  been  assumed  that  the  array  ins  been  aimed  in  the  direction  of  the 
signal  by  introducing  a  delay  of  n  seconds  per  hydrophone  where  i  is  the 
sampling  interval.  In  the  absence  of  a  signal  (the  hypothesis  H)  the  inputs 
to  the  M  hydrophones  are  all  uncorrelated  aid  arc  comprised  of  the  sum  of 
two  statistically  Independent  processes  n(t)  and  i(t)  with  powers  N  and  I 
respectively.  Under  the  alternative  (K)  ,  there  is  a  common  additive  signal 
s(t)  with  variance  S. 

The  expected  value  of  the  test  statistic  is  given  by 


E (S)  =  Tf  (M2-M) 
s 


jj^OOXj  (k) 


M  M 

h  Y  ^j(k>xr[k+m<r-l):  -  n  Y.  xj(k)xr[k4m(r“J)j 

r=l  r=l 


(A-2) 
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where  land  j  are  constrained  to  be  unequal.  Recognizing  that  only  those 


samples  at  the  same  instant  of  time  are  correlated  and  by  the  amount  2/tt 


sin  S/N+I+S,  this  mean  becomes, 


E(S)  =  Tf  (M2-M)  -  sin  1 

S  TT 


- § -  +  A  ) 

N+I+S  v  M  ..2' 
M 


(A-3) 


It  follows  that  the  shift  in  the  mean  of  the  test  statistic  becomes 
asymptotically  as  the  signal-to-noise  ratio  S/N  gets  small 

Ek(S)  -  EH(S)— T£s  <«2-«)  <1  -±)  jjfj- 


(A-4) 


Let  us  now  replace  the  quantity  (x^(k)  -  y^)  of  Eq .  (A-l)  by  the  symbol 
z^(k).  The  variance  of  the  test  statistic  under  the  hypothesis  can  be  written 


V“H(S)  *  Z_, 


Tf  Tf  M  M  M  M 

s  — ,  . — , _ 

,  '?  z  (k)z  (k)z  (p)z  (p) 

l—  i _ ,  : _ ,  i  J  qrnr 

k=l  p=l  i^j  q^n 


(A-5) 


The  z  terms  are  all  pair-wise  uncorrelated  unless  the  subscripts  differ  by 
the  right  amount  to  place  them  on  a  common  interference  front.  Therefore 

M.M 


Tf  M.M  Tf  /  :n 

-,s  - —  „  s 


Var^ (S)  « 


2Z  I  z2(k)  z2(k)+4  y  (Tfg-  m)  z±  (  k)  zi+^(k+mp)  (  k)  zj+p(  Wrap) 

k=l  i^j  p  =  1  i?fj  (A-6) 


Realizing  now  that  there  are  only  (M-p) (M-p-1)  terms  that  can  exist  in  the 
second  sura,  the  variance  becomes 

_  Tf  /  m  , 

,2 


VarH(S)  =  2Tfs  (M  -M)  z£  +  2  ^  (Tf^pm)  (M-p)  (M-p-1)  z±  (k)  (k-fm  p) 

p=i 

This  can  be  simplified  for  large  decision  times  by  recognizing  that 

2^jH-p)  (M-p-1)  =  (M2-M)  M 

where  M  =  (2M-4)/3.  Hence  for  large  decision  times 

-  —  2  2 


(A-7) 


Var..  (S)  c?  2Tf  (M  -M)  z.  (k)  +M  z,(k)  z _  (k-hnp ) 
H  s  i  i  i+p 


(A-8) 
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Recognizing  that 


(  1  ‘  if  i"J 


x1(k)xj  fk+n(j-l)]  =■  2 / 7T  8in_x  I/N+I',  if  j=i+p, 

0 •  otherwise 

\ 

it  is  straightforward  to  3how  that 


(A-10) 


2,.  ,  M-l 

Zi  00  “  IT 


1  -  —  sin 

7T 


-1  -i_l. 

N+I  I 


and 


Zi(k)z1+p(k4np)  =  - 

Combining  Eqs .  (A-10)  and  (A-ll) ,  it  is  seen  that 


2  -1  1  . 
—  sin  —  -  1 
it  N+l 


Var  (S)  =  2Tf  (M  -M) 
H  s 


2,-11  . 
-  sin  rr~~  -  1 
tt  N+I 


!  _  A  _  _i_l 

3M 


(A-ll) 


(A-12) 


The  output  signal-to-noise  ratio  of  the  detector  for  weak  signals  is  given  by 
the  ratic  of  the  shift  in  the  mean  to  the  square  root  of  the  variance  or 

2 

(A-13) 
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Appendix  B 

Let  pj,  denote  the  probability  that  a  sample  ir  contained  in  a  quantization 
interval  IK  given  that  the  interference  itself  is  contained  in  t'his  interval  •  i .  e . 

PK  “  Probj^n  I-  icl^/iel^  .  (B-l) 


If  the  interval  is  closed  and  of  width  A 
kA 


^(k-l)A  <  IK  <  kAj  , 


then 


PK  = 


J  |^kA-u j  -  ^(k-l)A-ujjf  i(u)  du  J^Fn(x)  -  Fn(x-A)j  f^kA-x)  dx 


_  (k-DA 


p(ieIK) 


f^(kA-x)  dx 
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If  pgc  is  the  probability  of  a  '  switch"  or  the  probability  that  a  sample 


lies  outside  this  interval  given  1  el,, ,  then 

In. 


p  00  -  1  -  pk 
sc 


a 

-'r\ 


(kA-x)dx  +  /  F  (x-A)  f.(A~x)dx 
n  I 


(B-3) 


P(1EIK) 


As  the  inteference-  to-noise  power  ratio  gets  large,  this  probability  of  a 
switch  becomes 


Psc<k> 


I/N  — £  oo 


I*  -  ■.»>) 


dx 


(B-A) 


Assuming  gaussian  processes,  this  can  be  put  into  the  form 


l(k  A)2  1  (k-l)A2  A//2tT 

i  IIP  exP  ~  2  t  exP  ~  2  t  r  r 

psc(k^ Hy - —  - - - — -1 —  *  /  erfc  xdx  <b-s> 


2p(iElK) 


0 


The  probability  of  a  switch  from  any  of  the  2 ( q  —  1)  closed  intervals  is  giv  i  by 

(B-6) 


(q— 1) 

"sc'1  -  2  £  k.c(k)'’(ltIK) 

k°l 


By  normalizing  the  quantization  intervals  by  the  standard  deviation  of  the 
noise ^ A  «=  m'^Nj  j  p  becomes 


Psc"^ 


n/ (2 

(q-1) 

{ 

nl 

j  erfc  xdx  < 

exp 

1  ,  2  2  N 
"  2  k  “  I 

+  exp 

-  \  (k-l)2m2  y 

J 

0 

k=0 

,  J 

1  In 

iT  Vi  11 

Invoking  symmetries,  thi9  can  be  written  as 


(B-7) 


P  ^ 
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7\|f  f1  /erfc  -d* 
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Ek  eXp 
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The  probability  of  a  switch  given  that  the  interference  lies  in  one  of 
the  two  open  intervals  is  given  by 


oo 

i 

/ 


psd(,>  -  ^ 


F^Cq-lM  -  u]f1(u),u  J[  1  -  Fn(x)j  fJx+Lq-DAjdx 


P  (i*-  I  ) 

q 


P(icIq} 


This  asymptotically  approaches 


PSO(^-;paT)J'  I1 
q 

which  becomes  for  gaussian  processes 


F  (x)  dx 
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(I 


(E 


_  1  .  .2  N  oo 

,  ,  ,  1  ,/n'  exP  -  2  (q_1)  T 

pso(q)‘~7  v  l/i  ^  *  /  erf c  xdx 
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2p ( i EI  ) 
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03 


0 


Recognizing  thatyr  j  er‘c  xdx  =  1,  the  probability  of  switching  from  either 
0 

of  the  open  intervals  becomes 
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